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ABSTRACT 

ZZ Ceti stars form the most numerous group of degenerate variable stars. They are 
otherwise normal DA (H-rich atmospheres) white dwarfs that exhibit pulsations. Here, 
we present an asteroseismological analysis for 44 bright ZZ Ceti stars based on a new 
set of fully evolutionary DA white dwarf models characterized by detailed chemical 
profiles from the centre to the surface. One of our targets is the archetypal ZZ Ceti 
star G117— B15A, for which we obtain an asteroseismological model with an effective 
temperature and a surface gravity in excellent agreement with the spectroscopy. The 
asteroseismological analysis of a set of 44 ZZ Ceti stars has the potential to characterize 
the global properties of the class, in particular the thicknesses of the hydrogen envelope 
and the stellar masses. Our results support the belief that white dwarfs in the solar 
neighbourhood harbor a broad range of hydrogen-layer thickness. 
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1 INTRODUCTION 

Pulsating DA (H-rich atmospheres) white dwarfs, commonly 
known as ZZ Ceti or DAV variable stars, comprise the most 
numerous class of compact pulsators. They are located in 
a narrow and probably pure instability strip with effective 
temperatures between 10 500 K and 12 500 K (e.g., Winget 
& Kepler 2008; Fontaine & Brassard 2008; Althaus et al. 
2010a). ZZ Ceti stars are characterized by multimode pho- 
tometric variations of up to 0.30 mag caused by nonradial 
g-mode pulsations of low degree {I ^ 2) and periods between 
70 and 1500 s. The driving mechanism thought to excite the 
pulsations is a sort of combination of the k — 7 mechanism 
acting in the hydrogen partial ionization zone (Dolez & Vau- 
clair 1981; Dziembowski & Koester 1981; Winget et al. 1982) 
and the "convective driving" mechanism proposed first by 
Brickhill (1991) and later re-examined by Goldreich & Wu 
(1999). The later mechanism is supposed to be dominant 
once a thick convection zone has developed at the stellar 
surface. 

White-dwarf asteroseismology fully exploits the com- 
parison between the observed pulsation periods in white 
dwarfs and the periods computed for appropriate theoreti- 
cal models. It allows us to infer details of the origin, internal 
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structure and evolution of white dwarfs (Winget & Kepler 
2008; Fontaine & Brassard 2008; Althaus et al. 2010a). In 
particular, constraints on the stellar mass, the thickness of 
the outer envelopes, the core chemical composition, weak 
magnetic fields and slow rotation rates can be inferred from 
the observed period patterns of ZZ Ceti stars. In addition, 
asteroseismology of ZZ Ceti stars is a valuable tool for study- 
ing axions (Isern et al. 1992; Corsico et al. 2001; Bischoff- 
Kim et al. 2008b; Isern et al. 2010, Corsico et al. 2011), crys- 
tallization (Montgomery & Winget 1999; Corsico et al. 2004, 
2005; Metcalfe et al. 2004; Kanaan et al. 2005), and impor- 
tant properties of the outer convection zones (Montgomery 
2005ab, 2007). Finally, the temporal changes in the observed 
stable periods allow the measurement of the white dwarf 
evolutionary timescale and the detection possible planets 
orbiting white dwarfs (MuUally et al. 2008). 

Among the numerous ZZ Ceti stars currently known 
(148 stars; Castanheira et al. 2010), in this paper we will 
analyze 44 bright ZZ Ceti stars which are listed in Table 1 
of Fontaine & Brassard (2008). We defer to a future work 
the study of the fainter ZZ Ceti stars discovered within the 
Sloan Digital Sky Survey (SDSS) (Mukadam et al. 2004; 
MuUally et al. 2005; Kepler et al. 2005b; Castanheira et al. 
2006, 2007, 2010). The first target star of our seismological 
survey is the most studied member of the class, the paradig- 
matic star G117— B15A. This star is an otherwise typical 
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DA white dwarf, the variabihty of which was discovered by 
McGraw & Robinson (1976) and, since then, it has been 
monitored continuously. The surface gravity, total mass, and 
effective temperature of this star have been the subject of 
numerous spectroscopic determinations. For instance, val- 
ues of logg = 7.97 ± 0.05, M* = 0.59 ± 0.03 M©, and 
Tcft = 11 630 ± 200 K, have been derived by Bergeron et al. 
(1995a, 2004) from optical spectra. Koester & Allard (2000) 
have reported somewhat lower values for the gravity and 
mass, logg = 7.86±0.14, M* = O.53±O.O7M0, and a higher 
effective temperature, T^s = 11 900 ± 140 K, from HST UV 
spectra. G117— B15A has oscillation periods 11 (amplitudes 
A) of 215.20 s (17.36 mma), 270.46 s (6.14 mma) and 304.05 
s (7.48 mma) (Kepler et al. 1982) that correspond to genuine 
pulsation modes. The star also shows the harmonic of the 
largest amplitude mode and two linear combinations. Kepler 
et al. (2005a) used the rate of change of the 215 s periodicity 
to show that the star has a C-0 core. The first detailed aster- 
oseismological study of this star was presented by Bradley 
(1998). This author obtained two different structures for the 
star according to the assignation of the radial order (fc) of the 
modes exhibited by the star. If the periods at 215, 271, and 
304 s are associated with k — 1,2,3, respectively, this au- 
thor obtained an asteroseismological model with a hydrogen 
envelope mass Mh/M* ~ 3 x 10~^. If, instead, the periods 
have k = 2,3, 4, the asteroseismological model was charac- 
terized by Mh/M, ~ 1.5 X 10~"'. Note that there are three 
orders of magnitude of difference in the mass of the H en- 
velope between the two possible (and nearly equally valid 
within their models) asteroseismological solutions. A sim- 
ilar degeneracy of seismological solutions for G117— B15A 
was also found by Benvenuto et al. (2002) on the basis of 
independent stellar and pulsation modeling. More recently, 
Castanheira & Kepler (2008) have found a seismological so- 
lution with Mh/M» ~ 10"'^ and k = 1,2,3 and another 
equally valid solution with Mh/M« ~ lO"'^ and k = 2,3,4. 
Finally, Bischoff-Kim et al. (2008a) also found two classes 
of solutions, one characterized by "thin" H envelopes, and 
other associated with "thick" H envelopes, although their 
"thick" envelope solutions (Mh/M» = 6 x 10~^) are consid- 
erably thinner than those of the previous works. 

Each of the mentioned asteroseismological studies con- 
stitutes a clear demonstration of the formidable capability 
of asteroseismology to shed light on the internal structure 
of DA white dwarfs. However, as important as they are, all 
of these studies are based on DA white-dwarf models that 
lack a fully consistent assessment of the internal chemical 
structure from the core to the outer layers. For instance, in 
the models of Bradley (1998), ahhough the C/He and He/H 
chemical interfaces are more realistic than previous stud- 
ies that used the trace element approximation (Tassoul et 
al. 1990), the core C-0 chemical profiles have a (unrealis- 
tic) ramp-like shape. In the case of Benvenuto et al. (2002) , 
the artificially-generated models are characterized by He/H 
chemical interfaces resulting from a time-dependent element 
diffusion treatment, and the C-O core chemical structure 
is extracted from the independent computations of Salaris 
et al. (1997). So, there is no consistent coupling between 
the chemical structure of the core and the chemical strati- 
fication of the envelope of the models. On the other hand, 
the recent works by Castanheira & Kepler (2008, 2009) are 
based on DA white-dwarf models similar to those of Bradley 



(1998), with a parametrization that mimics the results of 
time-dependent diffusion computations for the He/H chem- 
ical interfaces, but with a simplified treatment of the core 
chemical structure (50 % O and 50 % C). Finally, the study 
of Bischoff-Kim et al. (2008a) employs DA white-dwarf mod- 
els similar to those of Castanheira & Kepler (2008, 2009), 
but the envelope is stitched to a core that incorporates chem- 
ical profiles similar to those of Salaris et al. (1997). 

Needless to say, white-dwarf stellar models with consis- 
tent and detailed chemical profiles from the centre to the 
surface are needed to correctly assess the adiabatic pulsa- 
tion periods and also the mode-trapping properties of the 
DAVs, the crucial aspects of white-dwarf asteroseismology 
(Bradley 1996; Corsico et al. 2002). In this regard, Althaus 
et al. (2010b) (see also Renedo et al. 2010) have recently pre- 
sented the first complete set of DA white-dwarf models with 
consistent chemical profiles for both the core and the enve- 
lope for various stellar masses appropriate for detailed aster- 
oseismological fits of ZZ Ceti stars. These chemical profiles 
are computed from the full and complete evolution of the 
progenitor stars from the zero age main sequence, through 
the thermally-pulsing and mass-loss phases on the asymp- 
totic giant branch (AGB), and from time-dependent element 
diffusion predictions during the white-dwarf stage. 

In this paper, we carry out the first asteroseismologi- 
cal application of the DA white-dwarf models presented in 
Althaus et al. (2010b). Specifically, we perform a detailed 
asteroseismological study on 44 ZZ Ceti stars that includes 
the archetypal star G117— B15A, by using a grid of new evo- 
lutionary models characterized by consistent chemical pro- 
files and covering a wide range of stellar masses, thicknesses 
of the hydrogen envelope and effective temperatures. The 
asteroseismological analysis of such a large set of stars is a 
good starting point for ensemble asteroseismology of ZZ Ceti 
stars (see Castanheira & Kepler 2009). We also explore, in 
the frame of standard evolutionary calculations for the for- 
mation of DA white dwarfs, to what extent the mass of the 
He-rich envelope (Mhc) expected in DA white dwarfs de- 
pends on the details of prior evolution of progenitor stars. 
The paper is organized as follows. In Sect. (2] we provide a 
brief description of the evolutionary code, the input physics 
adopted in our calculations and the grid of models employed. 
There, we also explore the dependence of Mhc on the pro- 
genitor evolution. In Sect.O we describe our asteroseismo- 
logical procedure. In Sect. U we present our results, starting 
with the asteroseismological analysis for G117— B15A and 
a comparison with previous results (Sect. H?^ . and then by 
describing the results for the set of 44 stars (Sect. We 
conclude in Sect. [5] by summarizing our findings. 



2 NUMERICAL TOOLS AND MODELS 

2.1 Evolutionary code and input physics 

The present asteroseismological study is based on the full 
DA white dwarf evolutionary models of Althaus et al. 
(2010b) (see also Renedo et al. 2010) generated with the 
LPCODE evolutionary code. To our knowledge, these models 
are the first complete set of DA white dwarfs models char- 
acterized by consistent chemical profiles for both the core 
and envelope. This feature renders these models particularly 
suitable for asteroseismological studies of DA white dwarfs. 
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Here, we will briefly outline the most relevant charac- 
teristics of our evolutionary models of relevance for their 
pulsation properties. Further details can be found in Al- 
thaus et al. (2010b). In our computations, the ^^C(a,7)^®0 
reaction rate, of special relevance for the C-O stratification 
of the emerging white dwarf, is taken from Angulo et al. 
(1999). Thus, our white dwarf models are characterized by 
systematically lower central O abundances than the values 
predicted by Salaris et al. (1997), who used the larger rate 
of Caughlan et al. (1985). For example, for a ~ O.BIM© 
white dwarf, our computations give Xq ~ 0.73, about 4% 
lower than quoted by Salaris et al. (1997) (Xq ~ 0.76). Ex- 
tra mixing episodes during core He burning, of relevance for 
the final chemical stratification of white dwarfs (Straniero et 
al. 2003), was allowed to occur following the prescription of 
Herwig et al. (1997). Breathing pulses, which are convective 
runaways occurring towards the end of core helium burning, 
were suppressed. An important feature of our computations 
is that extra-mixing episodes were disregarded during the 
thermally-pulsing AGB phase, in line with theoretical and 
observational evidence (Herwig et el. 2007, Lugaro et al. 
2003, Salaris et al. 2009). This leads to the inhibition of the 
occurrence of the third dredge-up in low-mass stars, and 
consequently, to the gradual increase in the hydrogen-free 
core (HFC) mass as evolution proceeds during this phase. 
As a result, the initial-final mass relationship by the end of 
the thermally-pulsing AGB is markedly different from that 
resulting from considering the mass of the HFC right before 
the first thermal pulse. This issue is relevant for the C-O 
composition expected in a white dwarf. Depending on the 
white dwarf mass, the central oxygen abundance may be un- 
derestimated by about 15 % if it is assumed that the white 
dwarf mass is the HFC mass by the first thermal pulse (see 
Althaus et al. 2010b). 

We considered mass-loss episodes during the core he- 
lium burning stage and on the red giant branch following 
Schroder & Cuntz (2005) . During the AGB and thermally- 
pulsing AGB phases, we adopted the maximum mass loss 
rate between the prescription of Schroder & Cuntz (2005) 
and that of Vassihadis & Wood (1993). 

In our evolutionary computations, we have consid- 
ered the distinct physical processes that are responsible for 
changes in the chemical abundance distribution during white 
dwarf evolution. This is one of the most important improve- 
ments of our computations in comparison with previous as- 
teroseismological works on DA white dwarfs. In particular, 
element diffusion strongly modifies the chemical composi- 
tion profile throughout their outer layers in the course of 
evolution. As a result of diffusion processes, our sequences 
develop pure H envelopes and modifies the various inter- 
shells above the C-O core. We have considered gravitational 
settling as well as thermal and chemical diffusion — but not 
radiative levitation, which is relevant at high effective tem- 
peratures for determining the surface composition — of ^H, 
^He, *He, ^^C, "C, and ^''O (see Affhaus et al. 2003 for 
details). The standard mixing length theory for convection 
— with the free parameter a = 1.61 — has been adopted. 
Our treatment of time-dependent diffusion is based on the 
multicomponent gas treatment presented in Burgers (1969). 
In LPCODE, diffusion becomes operative once the wind limit 
is reached at high effective temperatures (Unglaub & Bues 
2000). In addition, abundance changes resulting from resid- 
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Figure 1. The grid of DA white dwarf evolutionary se- 
quences considered in this study represented in the plane M* — 
log{Mii/Mt). Each small circle corresponds to a sequence of DA 
white dwarf models with a given stellar mass and thickness of H 
envelope. The circles connected with a thick (orange) line corre- 
spond to the values of the maximum thickness of the H envelope as 
predicted by our evolutionary computations. For each sequence, 
we have pulsationally analysed about 200 stellar models covering 
the effective temperature range of 14 000 — 9 000 K. 



ual nuclear burning have been taken into account in our sim- 
ulations. Finally, we considered the chemical rehomogeniza- 
tion of the inner carbon-oxygen profile induced by Rayleigh- 
Taylor (RT) instabilities following Salaris et al. (1997). 

An important feature of our models is the dependence 
on the stellar mass of the outer layer chemical stratifica- 
tion expected in ZZ Ceti stars. Indeed, for the more massive 
models, diffusion strongly modifies the chemical abundance 
distribution, eroding the thick intershell region below the He 
buffer by the time evolution has reached the domain of the 
ZZ Ceti instability strip (see Althaus et al. 2010b). This is 
in contrast with the situation encountered in our less mas- 
sive models (M* < 0.63Mq), where the intershell region is 
not removed by diffusion. This is because element diffusion 
is less efficient in less massive models (with the subsequent 
longer diffusion timescale) and also because the intershell 
is thicker in these models. Regarding white dwarf astero- 
seismology, these are not minor issues, since the presence of 
a double-layered structure in the helium-rich layers is ex- 
pected to affect the theoretical g-mode period spectra of ZZ 
Ceti stars. It is clear that white dwarf evolution computed 
in a consistent way with element diffusion as considered in 
this study is required for precise asteroseismology. 

2.2 About the He content of a DA white dwarf 
star 

In this section we show that, in the frame of standard evolu- 
tionary computations for the formation of DA white dwarfs, 
the He content of these stars cannot be substantially smaller 
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Table 1. The values of the stellar mass of our set of DA white dwarf models (upper row) and the mass of H corresponding to the different 
envelope thicknesses considered for each stellar mass. The second row shows the maximum value of the thickness of the H envelope for 
each stellar mass according to our evolutionary computations. 
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Figure 2. Upper panel: change in the He content in the region 
limited by the boundaries of the HeFC and the HFC during the 
thermally-pulsing AGB phase. Lower panel: the temporal evolu- 
tion of surface luminosity and H- and He-burning luminosities in 
solar units for our initially 2MQ-star star during the thermally- 
pulsing AGB phase. 



than that predicted by our calculations. To do this, we com- 
pute the evolution of a 2A'/Q-star from the ZAMS until the 
thermally-pulsing phase on the AGB. The only way we en- 
visage in which the star may experience a substantial de- 
crease in its content of He is by undergoing a large num- 
ber of thermal pulses. In order for the model star to expe- 
rience the largest possible number of thermal pulses, and 
thus, the content of He decreases as much as possible, we 
switched off mass loss during this stage in our evolution- 
ary code. The results of this experiment are depicted in 
Fig. [2] in which we show the He content in the region lim- 
ited by the boundaries of the He-free core (HeFC) and the 
HFC in terms of time during the thermally pulsing phase 
(upper panel), and the surface luminosity and the H- and 
He-burning luminosities for each pulse in that phase (lower 
panel). We stopped the experiment when the object expe- 
rienced about 30 thermal pulses, which is enough for our 



purposes. We found that the He content of the object de- 
creased from Muc/Mq = 3.34 x lO"'^ (before the first ther- 
mal pulse) to MHe/M0 = 8.6 x 10"^ (before the thirtieth 
thermal pulse). Thus, the decrease (in solar masses) of the 
He content of the HFC is of a factor 3.89. However, it should 
be kept in mind that this reduction is due mainly to the in- 
crease of the mass of the future white dwarf, that grows 
from MfwD = O.523M0 to AffwD = 0.7114Mq between the 
thermal pulses 1 to 30. 

Our experiment shows that the He mass left in a DA 
white dwarf could be as much as a factor 3 — 4 lower than 
the values predicted by standard evolutionary computations, 
but not 2 or 3 orders of magnitude lower, which would be 
necessary for g-mode periods to be substantially afi'ected. 
We conclude that we can safety ignore the variation of Mhc 
in our asteroseismological analysis of ZZ Ceti stars. 

2.3 The model grid 

The DA white dwarf models employed in this study are the 
result of full evolutionary calculations of progenitor stars for 
solar-like metallicity [Z = 0.01). The complete evolution of 
eleven evolutionary sequences with initial stellar mass in the 
range 1 — SM© has been computed from the ZAMS through 
the thermally-pulsing and mass-loss phases on the AGB and 
finally to the domain of planetary nebulae. The values of 
the stellar mass of our set of models is shown in the upper 
row of Table [T] The range of stellar mass covered by our 
computations comfortably accounts for the stellar mass of 
most of the observed pulsating DA white dwarfs. 

Our asteroseismological approach basically consists in 
the employment of detailed white dwarf models character- 
ized by very accurate physical ingredients. These models are 
obtained by computing the complete evolution of the pro- 
genitor stars. We have applied successfully this approach to 
the hot DOVs or GW Vir stars (see Corsico et al. 2007a, 
2007b, 2008, 2009). Since the final chemical stratification of 
white dwarfs is fixed in prior stages of their evolution, the 
evolutionary history of progenitor stars is of utmost impor- 
tance in the context of white dwarf asteroseismology. Our 
asteroseismological approach, while being physically sound- 
ing, is by far much more computationally demanding than 
other approaches in which simplified models are used. As a 
result, our approach severely limits the exploration of the 
parameter space of the models. Indeed, for the case of DA 
white dwarfs, we have only two parameters which we are 
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able to vary in a consistent way: the stellar mass (M.) and 
the effective temperature (Tefi). Instead, the thickness of 
the H envelope (Mh), the content of He (A/hg), the shape of 
the C-O chemical structure at the core (including the pre- 
cise proportions of central O and C), and the thickness of 
the chemical transition regions are fixed by the evolutionary 
history of progenitor stars. Therefore, to push on the lim- 
its of our asteroseismological exploration, it would be de- 
sirable to change some additional parameters besides the 
stellar mass and effective temperature of our DA models. 
In this work, we have chosen to vary the thickness of the 
H envelope, because of the uncertainties in the mass loss 
rates. According to full evolutionary computations (Althaus 
et al. 2010b), the maximum H envelope mass expected in 
a white dwarf depends on the stellar mass and ranges from 
Mh/M. ~ 2.4 X 10"" (for M, = O.525M0) to 8.5 x lO"'' (for 
M* = O.878M0)(see the first row of Table [T]). Our decision 
for changing this parameter is due to several reasons: first, 
there are compelling theoretical reasons to believe that the 
H-content of DA white dwarfs might depend on the details 
of their previous evolution. On the contrary, the He content 
or the inner C-O chemical profiles are not expected to vary 
significantly due to the details of the previous evolutionary 
history (with the exception of a possible merger origin for 
the white dwarfs). Indeed, the total H content remaining in 
some DA white dwarfs could be several orders of magnitude 
lower than that predicted by our standard treatment of pro- 
genitor evolution. For instance, Althaus et al. (2005b) have 
found that Mh becomes considerably reduced if the progen- 
itor experiences a late thermal pulse episode (LTP) shortly 
after the departure from the thermally-pulsing AGB phase. 
In this sense, Tremblay & Bergeron (2008) show that the in- 
crease in the ratio of He-rich to H-rich white dwarfs can be 
understood on the basis that a fraction of DA white dwarfs 
above TcH ~ 10 000 K are characterized by a broad range of 
H-layer thickness. Second, the precise location of the He/H 
transition region (and the value of Mh) strongly affects the 
structure of the adiabatic period spectrum in a DA white 
dwarf (Bradley 1996). Finally, Mu is the structural param- 
eter that can be more easily modified in our models without 
removing relevant features predicted by the complete pro- 
genitor evolution. 

In order to get different thicknesses of the H envelope, 
we have followed a simple recipe. For each sequence char- 
acterized by a given stellar mass and a thick value of Mh, 
as predicted by the full computation of the pre-white dwarf 
evolution (second row of Table [T]), we have simply replaced 
^H by *He at the basis of the H envelope. This is done 
at very high effective temperatures (> 70 000 K), in such 
a way that the unphysical transitory effects associated to 
this procedure end much long before the models reach the 
stage of pulsating DA white dwarfs. After our ad hoc pro- 
cedure to change the thickness of the H envelope, we allow 
time-dependent element diffusion to operate while the mod- 
els cool down until they reach the effective temperatures 
characterizing the DAV instability strip. Diffusion leads to 
very smooth chemical profiles at the He/H chemical transi- 
tion regions. The resulting values of the H content for the 
different envelopes are shown in Table [1] and a graphical 
representation of the basic grid of models employed in this 
work is displayed in Fig. [1] In this figure, the canonical val- 
ues of AIh predicted by stellar evolution are connected with 



a thick (orange) line. Obviously, beyond the availability of 
the models of this coarse grid, we have the capability to 
generate additional DA white dwarf evolutionary sequences 
with arbitrary values of Mh for each stellar mass in order 
to refine the model grid. 



2.4 Pulsation computations 

We carried out the adiabatic pulsation computations re- 
quired by the present asteroseismological analysis by em- 
ploying the nonradial pulsation code described in Corsico & 
Althaus (2006). Briefly, the code, which is coupled to the 
LPCODE evolutionary code, is based on the general Newton- 
Raphson technique and solves the full fourth-order set of 
equations governing linear, adiabatic, nonradial stellar pul- 
sations following the dimensionless formulation of Dziem- 
bowski (1971). The prescription used to assess the run of 
the Brunt- Vaisala frequency (A'^) is the so-called "Ledoux 
modified" treatment (see Tassoul et al. 1990) appropriately 
generalized to include the effects of having three nuclear 
species varying in abundance. 

In the lower panel of Fig. [3] we show the spatial run of 
the logarithm of the squared Brunt- Vaisala frequency for 
models with M, = O.6O9M0 and different values of the 
thickness of the H envelope for Teft « 12 000 K. In the 
upper panel, we plot the internal chemical stratification of 
the models for the main nuclear species. The figure empha- 
sizes the role of the chemical interfaces on the shape of the 
Brunt- Vaisala frequency. In fact, each chemical transition 
region produces clear and distinctive features in A'^, which 
are eventually responsible for the mode-trapping properties 
of the models. In the core region, there are several peaks at 
— log(g) « 0.4 — 0.5 (where q = 1 — Mr/Mt^) resulting from 
steep variations in the inner C-O profile. The step shape of 
the C and O abundance distribution within the core, which 
is due to the occurrence of extra mixing episodes beyond the 
fully convective core during central helium burning, consti- 
tutes an important source of mode-trapping in the core re- 
gion — "core-trapped" modes (see Corsico & Althaus 2006). 
The extended bump in A''^ at — log(g) ~ 1 — 2 is another 
relevant source of mode-trapping. This feature is caused by 
the chemical transition of He, C and O resulting from nu- 
clear processing in prior AGB and thermally-pulsing AGB 
stages. It is worth noting that the shape of this transition is 
affected by diffusion processes which are operative at these 
evolutionary stages. Finally, there is the He/H transition re- 
gion, which is also another source of mode trapping, in this 
case associated with modes trapped in the outer H envelope. 

We have performed pulsation calculations on about 
(11 X 7 X 200) = 15 400 DA white dwarf models. In this 
account, we have considered the number of stellar mass val- 
ues (11), the number of thicknesses of the H envelope for 
each sequence (« 7), and the number of models (f» 200) 
with effective temperature in the interval 14 000 — 9 000 K, 
respectively. For each model, adiabatic pulsation g-modes 
with 1=1 and 1 = 2 and periods in the range 80 — 2000 s 
have been computed. This range of periods corresponds (on 
average) to 1 < A; < 50 for ^ = 1 and 1 < fc < 90 for £ = 2. 
So, more than ~ 2 x lO*" adiabatic pulsation periods have 
been computed in this work. 
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Figure 3. Upper panel: the internal chemical profiles of DA white dwarf models with M« = O.6O9M0, T^ff ~ 12 000 K, and different 
thicknesses of the H envelope. Only the main nuclear species are depicted. Lower panel: the run of the logarithm of the squared Brunt- 
Vaisala frequency for each model. Note the correspondence between the chemical transition regions (upper panel) and the resulting 
features in the shape of the Brunt- Vaisala frequency. For details, see the text. 



3 ASTEROSEISMOLOGICAL FITS 

We search for an asteroseismological model that best 
matches the pulsation periods of our target stars. To this 
end, we seek the model that minimizes a quality function 
defined simply as the average of the absolute differences be- 
tween theoretical and observed periods (e.g., Bradley 1998): 

JV 

cE. = $(M„ Mh, Teff) = ^ ^ - nt\ (1) 

!=1 

where N is the number of the observed periods in the star 
under study. We also have considered the quality function 
defined as (e.g., Corsico et al. 2009): 

JV 

= x'(M., Mh, Teff) ^^Y. ^^'^P*^' - 1^ (2) 

1=1 

Finally, we employ the following merit function (e.g., Cas- 
tanheira & Kepler 2008): 



H = H(M.,Mh,T.«) = ./ '"V"'T (3) 

i=i V 

where the amplitudes Ai are used as weights of each ob- 
served period. In this way, the period fit is more influenced 
by modes with large amplitudes than by the ones with low 
amplitudes. 

In the asteroseismological analysis of this work, we 
have employed the three quality functions $, , and H, 
deflned by Equations ([l}, ([2]), and ((3|, respectively. Since 
generally these functions lead to very similar results, we 
shall describe the quality of our period fits in terms of the 
function $ = $(M»,MH,Tcft) only. The effective temper- 
ature, the stellar mass and the mass of the H envelope 
of our DA white dwarf models are allowed to vary in the 
ranges: 14 000 ^ Tofi ^ 9 000 K, 0.525 ^ M, ^ O.877M0, 
—9.4 < log(A4"H/Mt) < —3.6, where the ranges of the values 
of Mh are dependent on M« (see Tabl^and Fig.[l|. For sim- 
plicity, the mass of He has been kept fixed at the value pre- 
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Figure 4. The location of the 44 ZZ Ceti stars analysed in this paper in the log g — T^g plane. The lines correspond to our set of DA 
white dwarf evolutionary tracks with thick (canonical) H envelope thickness. 



dieted by the evolutionary computations for each sequence. 
As we discussed in Sect. 12.21 the mass of the He content is not 
expected to be substantially smaller (say 100 — 1000 times) 
than predicted by our modeling. For this not too large uncer- 
tainty in the He content, only a weak dependence of the g- 
mode adiabatic pulsation periods on the value of Mhc is ex- 
pected (Bradley 1996), at variance with what happens with 
A/h. Finally, artificially changing the He mass of our mod- 
els would imply moving the triple transition C-O/He, which 
should introduce serious and undesirable artificial changes 
in the chemical structure of the models. The shape of the C- 
O chemical profile at the core and the central abundances of 
O and C have been also kept fixed according the predictions 
of the evolution during the central He burning stage of the 
progenitors. Finally, the thicknesses of the C-O/He and the 
He/H chemical transition regions have also been kept fixed 
at the values dictated by time-dependent element diffusion. 



4 STARS ANALYSED AND RESULTS 

We have carried out asteroseismological fits for a set of 44 
bright ZZ Ceti stars, the atmospheric parameters of which 
are shown in columns 2 and 3 of Table (2] In this Table, 
the stars have been sorted by decreasing Teg. The location 
of the studied stars in the loggr — T^s plane is displayed 
in Fig. [4] along with our evolutionary tracks. We defer to a 
future work the study of the fainter ZZ Ceti stars discovered 



within the SDSS. Most of these stars have been included in 
the study of Castanheira & Kepler (2009). 

In this section we present the results of our asteroseis- 
mological inferences. Because G117— B15A is the benchmark 
the ZZ Ceti class, we will devote the complete section 14.21 
to describe in detail the results of our asteroseismological 
analysis for this star, including a discussion of our findings, 
and defer the presentation of results for the whole sample 
of the analysed stars to the subsequent section. Before go- 
ing to the description of our asteroseismological results, we 
briefiy examine below the spectroscopic masses derived for 
the studied DAV stars and how the average value fits to the 
mean mass of DA white dwarfs reported by recent works. 



4.1 Spectroscopic masses 

The spectroscopic masses of the 44 ZZ Ceti stars studied in 
this work are shown in column 4 of Table (2] They have 
been derived simply interpolating from the tracks in the 
log g — Tefi diagram given the values of log g and T^g inferred 
from spectroscopic analysis. The mean value of the spectro- 
scopic masses for our sample of DAV stars is {M,}spac ~ 
0.630 ± O.O28M0. It is interesting to compare this value 
with the average mass of DA (pulsating and not pulsating) 
white dwarfs according to recent studies. Our value is some- 
what higher (~ 4%) than the value reported by Kepler et al. 
(2010) for DA white dwarfs on the basis of a large sample of 
1505 stars of the SDSS (DR4), (Af,)DA = 0.604 ± 0.003Mq , 
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Table 2. Atmospheric parameters and spectroscopic masses for 
the sample of ZZ Ceti stars analysed in this paper. 



Star 


Tefl [K] 


logs 


M,/Mq 


Ref. 


G226-29 


12 460 + 


200 


8.28 + 0.05 


0.771 +0.032 


3 


HS 1531+7436 


12 350 + 


181 


8.17 + 0.048 


0.704 + 0.029 


1 


G185-32 


12 130 ± 


200 


8.05 + 0.05 


0.634 + 0.028 


3 


L19-2 


12 100 ± 


200 


8.21 + 0.05 


0.726 + 0.033 


3 


G132-12 


12 080 + 


200 


7.94 + 0.05 


0.575 + 0.026 


4 


EC 11507-1519 


12 030 + 


200 


7.98 + 0.05 


0.596 + 0.026 


4 


PC 1541+650 


12 000 + 


70 


7.79 ± 0.04 


0.502 + 0.023* 


6 


R548 


11990 + 


200 


7.97 + 0.05 


0.590 + 0.026 


3 


CD 165 


11980 + 


200 


8.06 + 0.05 


0.639 + 0.029 


3 


GD 66 


11980 + 


200 


8.05 + 0.05 


0.634 + 0.028 


3 


G207-9 


11950 + 


200 


8.35 + 0.05 


0.812 + 0.033 


3 


EC 14012-1446 


11900 + 


200 


8.16 + 0.05 


0.696 + 0.031 


3 


KUV 11370+4222 


11890 + 


200 


8.06 + 0.05 


0.639 + 0.028 


3 


G238-53 


11890 + 


200 


7.91 + 0.05 


0.559 + 0.025 


3 


GD 99 


11820 + 


200 


8.08 + 0.05 


0.650 + 0.028 


3 


G29-38 


11820 + 


200 


8.14 + 0.05 


0.684 + 0.030 


3 


LP 133-144 


11800 + 


200 


7.87 + 0.05 


0.539 + 0.025 


3 


HS 1249+0426 


11 770 ± 


181 


7.92 ± 0.048 


0.564 + 0.024 


1 


MCT 2148-2911 


11 740 ± 


200 


7.82 + 0.05 


0.515 + 0.023* 


5 


GD 385 


11710 + 


200 


8.04 + 0.05 


0.627 + 0.028 


3 


GD 244 


11680 + 


200 


8.08 + 0.05 


0.650 + 0.028 


2 


HS 0507+0434B 


11630 + 


200 


8.17 + 0.05 


0.702 + 0.030 


3 


G117-B15A 


11630 + 


200 


7.97 + 0.05 


0.589 + 0.026 


3 


EC 23487-2424 


11520 + 


200 


8.10 + 0.05 


0.661 +0.028 


3 


MCT 0145-2211 


11500 + 


200 


8.14 + 0.05 


0.684 + 0.030 


3 


KUV 08368+4026 


11490 + 


200 


8.05 + 0.05 


0.633 + 0.028 


3 


PC 2303+243 


11480 + 


200 


8.09 + 0.05 


0.655 + 0.028 


3 


BPM 31594 


11 450 ± 


200 


8.11 + 0.05 


0.666 + 0.029 


3 


HLTau-76 


11450 + 


200 


7.89 ± 0.05 


0.548 + 0.025 


3 




11440 + 


200 


Q 1 7 + n 05 


n 709 + ft 030 


3 


HE 1429-037 


11434 + 


36 


7.82 + 0.02 


0.514 + 0.010* 


7 


G191-16 


11420 + 


200 


8.05 + 0.05 


0.632 + 0.028 


3 


HE 1258+0123 


11400 + 


200 


8.04 + 0.05 


0.627 + 0.029 


3 


G232-38 


11350 + 


200 


8.01 + 0.05 


0.610 + 0.027 


4 


KUV 02464+3239 


11290 + 


200 


8.08 + 0.05 


0.648 + 0.028 


2 


HS 1625+1231 


11270 + 


181 


8.06 ± 0.048 


0.638 + 0.027 


1 


BPM 30551 


11260 + 


200 


8.23 + 0.05 


0.737 + 0.032 


3 


HS 1824-6000 


11192 + 


181 


7.65 ± 0.048 


0.427 + 0.030* 


1 


G38-29 


11180 + 


200 


7.91 + 0.05 


0.557 + 0.025 


3 


GD 154 


11180 + 


200 


8.15 + 0.05 


0.689 + 0.029 


3 


R808 


11160 + 


200 


8.04 + 0.05 


0.626 + 0.028 


3 


BPM 24754 


11070 + 


200 


8.03 + 0.05 


0.620 + 0.028 


3 


G30-20 


11070 + 


200 


7.95 + 0.05 


0.578 + 0.026 


3 


PG 1149+058 


10 980 + 


181 


8.10 ± 0.048 


0.660 + 0.027 


1 



References: (i)Voss et al. (2006), (2) Fontaine et al. (2003). 
(3) Bergeron et al. (2004), t"*) Gianninas et al. (2006), 
(5)Gianninas et al. (2005), C^'Homeier et al. (1998),('') Silvotti et 
al. (2005). 

Note: the values of the stellar mass marked with * have been 
derived by extrapolation from our evolutionary model grid, and 
so, they are uncertain. 

and in agreement with the recent determination of Falcon 
et al. (2010), (M,)da = 0.647t;3 !^}^ Mq, obtained from the 
gravitational redshift determination of 449 DA white dwarfs, 
and that of Tremblay et al. (2011), {M,)da = O.eiSM©, us- 
ing 1089 DAs from DR4 of the SDSS. 



4.2 The archetypal ZZ Ceti star G117 B15A 

For this star, we initially computed the merit functions 
through our model grid by assuming that the harmonic de- 
gree of the three observed periods of G117— B15A is ^ = 1 
from the outset. Somewhat disappointing, we did not find 
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Figure 5. The quality function <I>(M*, Mjj, Toff) in terms of 
Tcff for the different values of the hydrogen thickness (shown 
with different colors and symbols) and a stellar mass of M, = 
O.593M0. The gray strip correspond to the effective temperature 
of G117— B15A according to spectroscopic analysis. Note in par- 
ticular the location of the best-fit solution (the minimum of <& at 
Tefi ~ 12 000 K), corresponding to the model 1 (^ = 1, fc = 2, 3, 4) 
in Table [3] 



any stellar model of the basic grid that matched simulta- 
neously the three observed periods. By closely examining 
our results, we discovered that a good period fit could be 
found by considering additional values of Mh near 
in the sequence with M» = 0.593Mq at approximately 
Teff = 12 000 K. Hence, we computed several additional se- 
quences with different values of Mn until a best-fit model 
with log(AfH/Af*) = —5.903 was found. The characteris- 
tics and periods of the best-fit model are shown in row 1 
(model 1) of Table [3] The period at 215 s has a radial or- 
der k = 2. Note that the fit to the main period is excel- 
lent (|A| = |n°''= - n"^| = 0.015 s), although the fits to 
the remainder two periods are not as good. The global fit, 
characterized by <E> = 1.729 s, is still very satisfactory. We 
repeated our computations but assuming 1=1 for the 215 
s mode at the outset, and allowing the other two periods to 
be associated with ^ = 1 or ^ = 2. We arrived at the same 
asteroseismological solution. 

In Fig. [S] we plot the function $(Af« , Mh, Toff ) in 
terms of the effective temperature for the different H enve- 
lope thicknesses corresponding to the sequence with M« = 
0.593AfQ. Clearly notorious is the existence of the best-fit 
solution at Toff - 12 000 K and Mh/M. « 1.25 x 10"''. 
Apart from the best-fit solution, there is another minimum 
at Tcff ~ 10380 K and Mn/M, ^ 1.17x 10"'', where $ ^ 4.5 
s. However, this solution must be discarded because its ef- 
fective temperature is too low as compared with the limits 
imposed by spectroscopy for G117— B15A. 

The uniqueness of the solution regarding the thick- 
ness of the H envelope is one of the main results of this 
work for G117— B15A. However, we warn that in this study 
we are matching 3 observed quantities (the pulsation pe- 
riods of G117— B15A) by varying just 3 structural quanti- 
ties {MtjMn, and Teg). So, it is not unconceivable that if 
we were varying an additional parameter (for instance Xq, 
Mhc, etc) of our models, we could found multiple asteroseis- 
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Table 3. Possible asteroseismological solutions for G117— B15A. Model 1 is the best-fit model corresponding to a family of solutions 
obtained by imposing that all of the observed periods correspond to £ = 1 modes. Models 2 to 4 result from the assumption that the 
observed periods are associated to £ = 1 or £ = 2 modes. 



Model 




M^/Mq 


log(MHo/A^.) 


log(MH/M,) 




TTttl 


£ 


k 


1 A 1 


<I> 




[K] 








[s] 


[s] 






[s] 


[s] 


1 


11986 


0.5932 


-1.62 


-5.90 


215.20 


215.215 


1 


2 


0.015 


1.729 












270.46 


273.437 


1 


3 


2.977 














304.05 


301.854 


1 


4 


2.196 




2 


12 450 


0.6090 


-1.61 


-4.45 


215.20 


214.947 


2 


6 


0.253 


0.177 












270.46 


270.268 


2 


8 


0.192 














304.05 


304.136 


2 


9 


0.086 




3 


12219 


0.6598 


-1.91 


-8.33 


215.20 


215.218 


2 


5 


0.018 


0.526 












270.46 


270.406 


1 


3 


0.054 














304.05 


305.557 


2 


8 


1.507 




4 


11735 


0.5930 


-1.62 


-9.33 


215.20 


214.422 


2 


4 


0.778 


0.735 












270.46 


271.682 


1 


2 


1.222 














304.05 


303.846 


2 


7 


0.204 





mological solutions due to the ambiguity introduced by the 
new parameter to be adjusted. Our models do not have an 
extra fit parameter. 

We also carried out additional period fits in which the 
value of £ for each of the theoretical periods is not fixed but 
instead is obtained as an output of our period fit procedure, 
although the allowed values are just 1=1 and £ — 2. The 
results are displayed in rows 2 to 4 of Table O For these 
models, the period fits are excellent. In particular, the peri- 
ods of model 2 match the observed periods with an average 
difference of ~ 0.18 s. One of the reasons is that are more 
£ = 2 modes per period interval. However, for the three 
models, the main periodicity of G117— B15A at 215 s is as- 
sociated to a ^ = 2 mode. This is in strong contradiction 
with the results of Robinson et al. (1995), who identify the 
215 s period with a, £ = 1 mode by means of time-resolved 
ultraviolet spectroscopy. This result is consistent with the 
further analysis of Kotak et al. (2004). Thus, as tempting 
as these solutions seem, they must all be discarded from our 
analysis. 



4-2.1 Estimation of the internal uncertainties 

We have assessed the uncertainties in the stellar mass (aM,), 
the thickness of the H envelope (ca/h) the effective tem- 
perature (o"T,jf) of the best-fit model by employing the ex- 
pression (Zhang et al. 1986; Castanheira & Kepler 2008): 



where So = <E>(M°, M^, T^g) is the minimum of $ which 
is reached at (M'^ , M^,T^g) corresponding to the best-fit 
model, and S is the value of $ when we change the parame- 
ter i (in this case, M, , Mh, or Toh) by an amount di, keeping 
fixed the other parameters. The quantity di can be evalu- 
ated as the minimum step in the grid of the parameter i. 
We obtain the following uncertainties, which are the inter- 
nal errors of our asteroseismic procedure: aM, ~ 0.007Mq, 
o-Mh ~ 0.7 X lO'^M*, and crT^f ~ 200 K. The uncertainties 



in the other quantities {L,,R,,g, etc) are derived from the 
uncertainties in Af, and T^b. 

In Table |4l we compare the main characteristics 
of our best-fit model with the observed properties of 
G117— B15A. In particular, we include the surface parame- 
ters of G117— B15A taken from several spectroscopic studies. 
We include also the spectroscopic mass computed by inter- 
polating from our evolutionary tracks. Note the agreement 
between the effective temperature and gravity of our aster- 
oseismological model and the values derived by Koester & 
AUard (2000) and Koester & Holberg (2001). The total mass 
of our model, however, is 3 — 11 % higher than the values de- 
rived in those studies. Our model is ~ 350 K hotter than the 
spectroscopic temperature of Bergeron et al. (1995a, 2004), 
and about 400 K cooler than the value derived by Robin- 
son et al. (1995), but the surface gravity and mass are in 
excellent agreement with the values quoted in both studies. 

4-2.2 Asteroseismological distance 

Since we have the luminosity of the best-fit model, we 
can estimate the asteroseismological distance and paral- 
lax of G117— B15A by means of the relation logci[pc] = 
1 (mv - My + 5) , where Mv = Afboi - BC- The bolometric 
magnitude, Afboi, can be computed as Afboi = Afboi(©) — 
2.5 log(L*/L0), being the bolometric magnitude of the Sun 
Mboi(0) = +4.75 (Allen 1973). By using BC = -0.611 
(Bergeron et al. 1995a) and mv = 15.50 (Bergeron et al. 
1995b), we obtain a distance d — 60.3 ± 2.5 pc, and a par- 
allax TT — 16.6 ± 0.8 mas, in excellent agreement with the 
inference of Bradley (1998) (tt = 16.5 mas). The distance 
estimated from optical, lUE, and HST spectra is 58 ± 2 pc, 
59 ± 5 pc, and 67 ± 9 pc, respectively. Holberg et al. (2008) 
derive a distance of 57.68 ±0.60 pc. Our seismological paral- 
lax is larger than the trigonometric value extracted from the 
Yale Parallax Catalog (van Altena et al. 1994) of 10.5 ± 4.2 
mas. In order for the asteroseismological parallax to be com- 
patible with the trigonometric one, the mass of the astero- 
seismological model should be as low as ~ O.35Af0! We can 
safely discard a low stellar mass for G117— B15A from spec- 
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Figure 6. The periods of the modes with i = 1 and k = 1, 2, 3 and 4 in terms of the thickness of the H envelope corresponding to 
Kh = O.525Af0 (left panel), Af* = O.593M0 (middle panel), and M* = O.877M0 (right panel). The observed periods of G117-B15A 
are shown with thin horizontal dashed lines. The dashed (orange) line in the middle panel indicates the match between the theoretical 
periods of the asteroseismological model and the periods observed in G117— B15A. 



troscopy. Then, we conclude that the trigonometric parallax 
must be more uncertain than quoted, and that the astero- 
seismological parallax is robust. 

4-2.3 Discussion 

All the previous asteroseismological studies on G117— B15A 
(Bradley 1998, Benvenuto et al. 2002, Castanheira & Ke- 
pler 2008, Bischoff-Kim et al. 2008a) report an ambiguity 
of the solutions regarding the thickness of the H envelope 
of this star. In those studies, a family of thin envelope so- 
lutions is obtained for a identification k — 1,2,3, whereas 
a second family of solutions of thick envelopes is derived 
if A: = 2, 3, 4. In contrast, our asteroseismological analysis 
strongly points to a single solution regarding the thickness 
of the H envelope, with a value of log(MH/M«) ~ —5.9, 
that corresponds to the identification k = 2,3, 4, which 
was associated to thick H envelopes in the previous stud- 
ies. The degeneracy of solutions is solved for the first time 
by our results. The reason for the uniqueness of the solu- 
tion in our computations is that, regardless of the value 
of the stellar mass, temperature, or thickness of the H en- 
velope, it is impossible to find a model whose mode with 
k = 1 has a period close to 215 s. This preclude us to find 
any possible asteroseismological model with the identifica- 
tion k = 1,2, 3. This is shown in Fig. [6l where we plot the 
periods in terms of the thickness of the H envelope for mod- 
els with Teff ~ 12 000 K and A/. = 0.525Mq (left panel), 
M. = O.593M0 (middle panel), and M. = 0.877 Mq (right 
panel). In the middle panel, the best-fit model is indicated 
with a dashed (brown) line. Clearly, the period for the mode 
with = 1 in our models is always very short in comparison 
with the shortest period shown by G117— B15A. We men- 
tion that we have also computed an additional evolutionary 
sequence with the same characteristics as our best-fit model, 
but with a thinner H envelope than that considered in Table 
[U (log(MH/M*) < -9.33). Even in this case, the pulsation 
periods exhibit the same trend shown in the middle panel 
Fig. [6l with the period of the k = 1 mode markedly de- 



parted from 215 s, and as a result, we are not able to find a 
thin-envelope solution. 

It should be kept in mind, however, that we could run 
into multiple solutions for G117— B15A if we were to vary an 
additional parameter of our models. For instance, it could 
be possible that if we were freely changing the He content 
(Mhb) of our models, for instance by adopting a He layer 
mass two order of magnitude thinner, then the k = 1 period 
could became close to 215 s, and so, we could recover the 
two families of thin and thick He envelope solutions found 
in the previous studies. However, as discussed in Sect. 12.21 
such low Mhc values are difficult to conceive from stellar 
evolution calculations. 

A distinctive feature shown in Fig. [6] is the presence of 
a behavior reminiscent to the well known "avoided crossing" 
(see also Fig. 3 of Castanheira & Kepler 2008). When a pair 
of modes experiences avoided crossing, the modes exchange 
their intrinsic properties (see Aizenman et al. 1977). In our 
models, avoided crossing is produced when we vary the 
thickness of the H envelope. As a result, for certain values of 
Mh, the period spacing turns out be very short. This effect 
is more notorious for low radial order modes. For instance, 
for the sequence with A/, — O.525A/0, the period spacing 
between the modes with k = 1 and fc = 2 is of only « 8 
s if log(AfH/Af*) ~ —4.8! Something similar is seen for the 
sequences with M, = O.593M0 and M, = O.877Af0, with H 
envelopes of log(AfH/M*) ~ -5.3 and log{M-a/Mt) ~ -7.4, 
respectively. We note also that avoided crossing is present in 
our models, but to a less extent, when we vary the effective 
temperature. 

In Figs. [7] and [S] we display the thin and thick envelope 
solutions with empty and filled symbols, respectively, for 
the asteroseismological solutions found in previous works. 
The location of our asteroseismological model is depicted 
with a (magenta) star symbol. Note that Bischoff-Kim et al. 
(2008a) found several equally valid asteroseismological mod- 
els with thin and thick H envelopes for G117-B15A. Fig. [7] 
shows a clear correlation between the mass and the effec- 
tive temperature of the solutions in the studies of Bradley 
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Table 4. Characteristics of G117— B15A and of our seismological model. The quoted uncertainties in the seismological model are the 
internal errors of our period-fit procedure. The progenitor star of the asteroseismological model star has a stellar mass of Af» = 1.75M0 
at the ZAMS. 



Quantity 


Robinson et al. 


Koester & AUard 


Koester & Holberg 


Bergeron et al. 


Our seismological 




(1995) 


(2000) 


(2001) 


(1995a, 2004) 


model 


Teff [K] 


12 375 ± 125 


11 900 ± 140 


12 010 ± 180 


11 630 ± 200 


11985 ± 200 


M,/M@ 


0.591 ± 0.031 


0.534 ± 0.072 


0.575 ± 0.092 


0.589 ± 0.026 


0.593 ±0.007 


logs 


7.97 ±0.06 


7.86 ±0.14 


7.94 ±0.17 


7.97 ±0.05 


8.00 ±0.09 


l0g(fl,/i?Q) 










-1.882 ± 0.029 


log(L./L0) 










-2.497 ± 0.030 












2.39 X 10-2 












(1.25 ±0.7) X 10-6 


Xc,Xo (center) 










0.28,0.70 



Note 1: the values of the spectroscopic mass quoted in columns 2 to 5 have been computed by interpolating from our set of 
evolutionary tracks (see Fig.[4]| using the corresponding values of logg and T^ff. 
Note 2: Robinson et al. (1995) use MLT/« = 1 model atmospheres, while the other use MLT/a = 0.6, hence they obtain lower Tegs. 



0.85 

0.8 
0.75 

0.7 
0.65 

0.6 
0.55 

0.5 
0.45 

0,4 



hollow: k= 1, 2, 3 ("thin") 
filled: k= 2, 3, 4 ("thick") 



Robinson et al. (1995) — 
Koester & Holberg (2001 ) 
Koester & Allard (2000) 
Bergeron et al. (1995. 2004) ~ 
Bradley (1998) 

Benvenuto et al. (2002) _ 
Castanheira & Kepler. (2008) 
Bisch^ff-Kim et al. (2008a) 
This work — 




_L 



_L 



_L 



_L 



10500 11000 11500 12000 12500 13000 13500 
TettIK] 

Figure 7. The location of the asteroseismological models for 
G117— B15A in the plane T^ff — M, according to the studies 
carried out up to date and according to this work, as indicated 
with different symbols. Empty symbols correspond to solutions 
for which the radial order identification of the observed periods 
is k = 1,2,3 (thin H envelopes), and filled symbols are associ- 
ated to solutions for which the observed periods have fc = 2, 3, 4 
(thick H envelopes). Also included is the location of G117— B15A 
according to several spectroscopic studies (filled squares). Note 
that Robinson et al. (1995) use MLT/a = 1 model atmospheres, 
while the other use MLT/o = 0.6, hence they obtain lower Tcfi 
values. 



(1998), BischoflF-Kim et al. (2008a) and Castanheira & Ke- 
pler (2008): cooler solutions have larger masses. The oppo- 
site trend is exhibited by the two solutions of Benvenuto et 
al. (2002). All the solutions of BischofT-Kim et al. (2008a) 
and two solutions of Castanheira & Kepler (2008) are sub- 
stantially more massive than the best-fit models of Bradley 
(1998), Benvenuto et al. (2002), and our own asteroseismo- 
logical model, and also than the estimations of the stellar 
mass of G117— B15A from spectroscopic studies. 



-10 



-9 



E -6 



1 I , I I I r 
hollow: k= 1, 2, 3 ("thin") 
filled: k= 2, 3, 4 ("thick") 
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Benvenuto et al. (2002) 
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This work _ 
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_L 



_L 



_L 



_L 



_L 
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0.65 

M./M 



0.7 0.75 0.8 0.85 



Figure 8. Similar to Fig.H but in the plane M, - log(AfH/A^'»)- 
The hollow green rectangle corresponds to the thin solutions of 
Bischoff-Kim et al. (2008a) corresponding to an identification k = 
1, 2, 3. In the interests of comparison, the location of the post-Late 
Thermal Pulse remnant with a thin H envelope of the scenario of 
Althaus et al. (2005b) is also shown. 



The existence of two separate families of solutions re- 
garding the thickness of the H envelope, as predicted by pre- 
vious studies, is clearly emphasized in Fig. [8] Here, it is no- 
table a correlation between A/h and M, , according to which 
the more massive asteroseismological models have thinner 
envelopes. Curiously, this trend is in line with the predic- 
tions of the canonical evolutionary computations (Althaus 
et al. 2010b), which are shown with a thick (orange) line 
that connects the maximum values of Mh for different stellar 
masses. The figure also shows a notable agreement between 
the thin and thick solutions of Bradley (1998), Benvenuto 
et al. (2002), and Castanheira & Kepler (2008). In contrast. 
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the solutions of Bischoff-Kim et al. (2008a) (both thin and 
thick) appear shifted toward smaller values of the H envelope 
thickness. In this context, our single seismological solution 
seems to be more nearly compatible with the family of thin 
H envelopes (although with a mode identification typical of 
thick-envelope solutions) than with the group of thick en- 
velopes. 

Although our new value for the H envelope thickness of 
G117— B15A is significantly lower than the canonical value 
predicted by stellar evolution (roughly 2 orders of magni- 
tude thinner), it is in perfect agreement with what could be 
expected from the LTP scenario. In this scenario, also called 
AGE final thermal pulse or AFTP scenario, a final helium 
shell flash is experienced by a star shortly after the depar- 
ture from the AGE (Blocker 2001). During a LTP, not all 
the hydrogen is burnt, in contrast to post-AGB stars that 
experience a very late thermal pulse (a born-again episode), 
but part is diluted by surface convection and mixed inwards 
with the underlying intershell region formerly enriched in 
helium, carbon and oxygen. Althaus et al. (2005b) have ex- 
plored the possibility that an initially 2.7Mq star experienc- 
ing a LTP shortly after the departure from the thermally 
pulsing AGE could reach the final cooling branch with a 
H envelope substantially smaller than predicted by stan- 
dard stellar evolution. They found that most of the origi- 
nal H-rich material of the post-AGB remnant is burnt after 
the post-LTP evolution, when the star returns to the high 
Teff regime for the second time, resulting in a white dwarf 
remnant of M» — 0.5885Mq with a value of the H enve- 
lope thickness of A/h = 1-7 x 10~®M,. Very interestingly, 
our best-fit model for G117— B15A and the DA white dwarf 
model resulting from the scenario proposed by Althaus et al. 
(2005b) are located roughly at the same place in the plane 
M* - log(MH/M*) (see Fig. UlE]. Therefore, our study rein- 
forces the validity of the results of Althaus et al. (2005b) 
about the existence of DA white dwarfs with H envelopes 
substantially thinner than the canonical value, and suggests 
that G117— B15A could be the descendant of a progenitor 
star that experienced a LTP episode before reaching the final 
cooling branch. 

4.3 The set of 44 ZZ Ceti stars 

Here, we present the asteroseismological analysis for the 44 
ZZ Ceti stars listed in Table [2] G117-E15A included. In 
the second and third columns of Table 5 we show the ob- 
served periods and amplitudes, respectively. These values 
are extracted from the works of Castanheira & Kepler (2008, 
2009), unless indicated otherwise. The fourth column of Ta- 
ble 5 shows the theoretical periods of the adopted asteroseis- 
mological model for each star, whereas the fifth and sixth 
columns include the i- and fc-identification, respectively, of 
each pulsation mode. The seventh column shows the ab- 
solute difference between observed and theoretical periods. 



We warn, however, that both models have diflfcrent internal 
chemical structure, in particular due to the presence of extra 
chemical structure and appreciable amounts of at the base of 
the He buffer in the post-LTP DA white dwarf model of Althaus 
et al. (2005b) (see the panel D of Figure 2 of Miller Bertolami et 
al. 2005). 



and the eighth column indicates the value of the quality 
function defined by Eq. IT}. 

Below, we describe the general criteria adopted to 
choose the asteroseismological model for each star. 



4-. 3.1 Criteria used in the fits 

Usually, when performing period-to-period fits to ZZ Ceti 
stars, we found multiple seismological solutions, that is, 
many stellar models that nearly reproduce the periods ob- 
served in a given DAV star. So, in order to isolate a single 
asteroseismological model among the several possible and 
equally valid ones, we must apply some criteria: 

• First, we looked for the models associated to the lowest 
value of the quality functions, thus ensuring that the ob- 
served periods are closely matched by the theoretical ones. 

• When possible, we used the external ^-identifications 
of the observed periods according to studies that employ 
the high-speed photometry method (see, e.g., Robinson et 
al. 1995), the time-resolved ultraviolet spectroscopy method 
(see, e.g., Kepler et al. 2000) or the time-resolved optical 
spectroscopy approach (see, e.g., Clemens et al. 2000). 

• When several families of solutions were found, we 
elected the models with values of Tch and logy as close as 
possible to the spectroscopic ones. In this way, we guaran- 
tee that the surface parameters of the asteroseismological 
solutions are not in conflict with observations. 

• Among possible asteroseismological solutions with simi- 
lar values of the quality function, we prioritized the solutions 
that fit the largest amplitude modes with theoretical modes 
having 1 = 1. This is because the well-known property that 
£ = 1 modes exhibit substantially larger amplitudes than 
£ — 2 ones, because geometric cancellation effects become 
increasingly severe as £ increases (Dziembowski 1977). 

• In the cases in which several modes had similar ampli- 
tudes in the power spectrum, we gave more weight to stellar 
models that fit those periods with theoretical periods hav- 
ing the same £ value. In this way, we are assuming that two 
eigenmodes with different values of the harmonic degree £ 
usually should not have similar amplitudes. 

• For a given star showing a large number of modes, we 
favored the seismological solutions that fit to observed peri- 
ods with a larger number of £ = 1 than £ — 2 modes. This is 
because there is more chance to observe £ — 1 modes than 
£ = 2 modes. 

• In the opposite case, for stars exhibiting just a single 
period, we employed only the set oi £ = 1 periods to per- 
form the period fit. Then, we chosen the asteroseismologi- 
cal model by searching for that model having the minimum 
value of the quality function, and we restricted the solu- 
tions by using the spectroscopic constraints {Tcb and logg), 
if necessary. 



4-3.2 Some particular cases 

Because the large number of ZZ Ceti stars seismologically 
analysed in this work, it would be unpractical and tedious to 
describe in detail the procedure we followed to arrive at the 
asteroseismological model for each star, as we already did 
for the particular case of G117— B15A. Instead, we briefiy 
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Table 5. Periods observed in the sample of 44 bright ZZ Ceti stars studied in this work and the corresponding theoretical periods and 
{£, fc)-identification of our asteroseismological models. 



Star 


obs 

n° = 


A 


TTth 


I 


k 








[s] 


[mma] 


[s] 










HS 1531+7436 


112.50 




112.499 


1 


1 


0.001 


0.001 


r^Tt 9/1/1 (1) 


zuz.yo 


/I ri/i 


lyo.y / o 


Z 


O 


7 nnv 
/ .UU r 


Z. 100 




256.56 


12.31 


257.215 


1 


3 


0.665 






294.60 


4.85 


296.820 


2 


9 


2.220 






307.13 


20.18 


306.283 


1 


5 


0.847 






QOf^ OS 
yuu.uo 


1 79 




\ 


1 Q 


u.uou 




G226-29 


109.28 




109.246 


1 


1 


0.032 


0.032 




Qf^c; fin 


0/1 n 

Z4.U 


ooO. r O i 


1 
i 







n VT'Q 




AA(K on 

440. ZU 


1 Q o 

io.y 


AACK A OO 
440. 4zy 


1 
i 


Q 
O 


n 090 
u.zzy 






555.30 


16.6 


556.767 


1 


11 


1.468 






743.40 


7.6 


742.920 


1 


16 


0.679 




I'D IQ"? 1/1/1 


ono on 


1 n 
iU 


01 1 O/IV 
Zl 1.Z4 / 


1 
1 


Z 


O C\A1 
Z.U4 ( 


i.ZoD 




Qn*^ vn 


o.o 


OU4.oy4 




z 


Q 


±.0U0 






327.30 


4.0 


327.716 


2 


9 


0.416 




EC 11507-1519 


191.70 


3.59 


191.964 


1 


2 


0.264 


0.231 




249.60 


7.70 


249.798 


1 


4 


0.198 






T 1 Q cn 


O A 
Z,4 




z 


Z 


n /I fiv 


1 90/1 
i.ZZ4 




1 1 a vn 


1 O 

i .z 


T 1 /I /loc: 
1 i4.4yo 


1 
1 


1 


A 9nc: 






1 A'^ fin 

140.DU 


U.D 


1 070 




z 


Q 


n 1 9S 

U. ±ZO 






192.60 


6.5 


192.561 


1 


2 


0.039 






350.10 


1.1 


351.359 


1 


6 


1.259 




^JJ DO ^ ' 




A 01 

4.Zi 


1 OQ 1 f\A 

lyo. iU4 


z 


4 


U.40U 






255.87 


3.43 


256.137 


2 


6 


0.270 






271.71 


16.70 


271.804 


1 


3 


0.089 






302.77 


11.29 


'^nn 1 09 

OUU. 


]^ 


4 


2.663 




G132-12 


212.70 


4.3 


212.703 


1 


2 


0.003 


0.003 


G207-9 


259.10 


17.3 


258.853 


1 


4 


0.247 


0.767 




292.00 


49.0 


290.379 


2 


10 


1.621 






318.00 


64.0 


318.257 


1 


5 


0.257 






557.40 


63.4 


556.204 


1 


12 


1.376 






740.40 


46.4 


741.034 


1 


17 


0.334 




G117-B15A 


215.20 


17.36 


215.215 


1 


2 


0.015 


1.729 




270.46 


6.14 


273.437 


1 


3 


2.977 






304.05 


7.48 


301.854 


1 


4 


2.196 




MOT 2148-2911 


260.80 


12.6 


260.798 


1 


4 


0.002 


0.002 



summarise below a few details related to the selection pro- 
cess of the best-fit model for some cases of interest. The 
structural parameters of the asteroseismological models for 
the complete set of ZZ Ceti stars analysed in this study are 
shown in Table [6l 

- G226-29. G226-29 also exhibits a single mode with 
a short period. Fortunately, there exist a robust constraint 
on its ^-identification. In fact, Kepler et al. (2005b) found 
that the mode is actually a triplet ((. = 1) with the central 
component at a period of 109.278 s. The solution in this 
case corresponds to a rather massive model with a thick H 
envelope (M. - O.77M0, Mh = 2.02x lO'^^M.), in line with 
the spectroscopic observations. 



- HS 1531-1-7436. This star exhibits a single mode with 
a very short period (for ZZ Ceti standards) at 112.5 s. Un- 
fortunately, the presence of just one period turns very dif- 
ficult any attempt of asteroseismology on this star, and we 
are forced to make a somewhat arbitrary assumption. If we 
assume that this mode corresponds to a (I, k) = (1, 1) iden- 
tification, then the stellar mass of the seismological model 
must be larger than 0.705Mq. In the first attempt to fit 
its periods, we obtained massive solutions (M« ~ 0.77Mq), 
but at effective temperatures excessively low (~ 10 800 K). 
These solutions are characterized by thick H envelopes. Since 
we have just a single observed period, it is possible to find a 
model with the appropriate H envelope thickness as to allow 
to fit the period at an effective temperature in close agree- 
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Table 5 - continued 



Star 


j^obs 




nth 




k 


lAI 


3> 




fsl 

PJ 


[mrrici] 












G38-29 


413.307 


3.07 


413.985 


2 


16 


0.678 


1.515 




432.354 


3.57 


434.227 


2 


17 


1.873 






546.960 


6.97 


545.442 


2 


22 


1.519 






705.970 


18.44 


707.049 


1 


16 


1.079 






840.390 


5.19 


839.307 


1 


20 


1.083 






899.971 


10.59 


896.903 


2 


38 


3.068 






922.567 


5.94 


921.066 


2 


39 


1.591 






945.448 


12.34 


946.328 


2 


40 


0.880 






962.007 


8.09 


962.277 


1 


23 


0.270 






963.593 


4.58 


962.277 


1 


23 


1.316 






yoy . / iy 


1 n fi/i 
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z 
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7.14 


1003.878 


1 


24 
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1016.15 


5.79 


1014.220 


2 


43 


1.930 






1081.82 


5.04 


1082.720 


1 


26 


0.900 




PG 1541+650 


689.00 




688.891 


1 


11 


0.109 


0.270 




(0 / .UU 




/0/.U4( 


i 


1 o 


U.U4 1 






564.00 




563.346 


2 


16 


0.654 




G191-16 W 


510.00 




509.983 


1 


9 


0.017 


0.931 




600.00 




598.812 


1 


11 


1.188 






710.00 




712.027 


1 


14 


2.027 






893.00 




893.495 


1 


18 


0.495 




G185— 32 


215.74 


1.93 


215.739 


1 


2 


0.001 


1.691 




266.17 


0.46 


269.253 


2 


7 


3.083 






300.60 


1.04 


298.724 


2 


8 


1.876 






370.21 


1.62 


367.694 


1 


5 


2.516 






651.70 


0.67 


652.677 


1 


12 


0.978 




EC 14012-1446 


398.90 


12.1 


403.823 


1 


7 


4.923 


2.541 




530.10 


16.7 


524.782 


1 


10 


5.318 






610.40 


54.3 


613.677 


1 


12 


3.277 






678.60 


7.6 


675.620 


1 


14 


2.980 






722.90 


22.9 


721.733 


1 


15 


1.167 






769.10 


51.7 


769.121 


1 


16 


0.042 






882.70 


2.9 


883.878 


2 


34 


1.178 






937.20 


11.0 


934.485 


2 


36 


2.715 






1217.40 


7.5 


1216.141 


1 


27 


1.259 




EC23487-2424 


804.50 


19.3 


806.160 


1 


19 


1.660 


2.297 




868.20 


12.8 


863.294 


1 


21 


4.906 






992.70 


24.4 


992.375 


1 


24 


0.325 




GD 165 


114.30 




114.278 


2 


2 


0.022 


0.889 




120.36 




119.195 


1 


1 


0.445 






192.68 




192.102 


1 


2 


0.578 






249.90 




252.412 


1 


3 


2.512 




R548 


187.28 


0.9 


187.597 


1 


1 


0.308 


2.516 




212.95 


5.4 


213.401 


1 


2 


0.451 






274.51 


3.5 


242.263 


1 


3 


2.249 






318.07 


1.1 


311.361 


2 


8 


6.709 






333.64 


1.3 


336.504 


2 


9 


2.864 




HE 1258+0123 


439.20 


9.8 


446.066 


2 


14 


6.867 


2.099 




528.50 


9.3 


527.704 


1 


9 


0.796 






628.00 


15.2 


627.326 


2 


21 


0.679 






744.60 


22.9 


744.780 


1 


14 


0.180 






881.50 


17.6 


892.728 


1 


17 


1.228 






1092.10 


14.1 


1094.947 


1 


22 


2.847 





GD 154 



402.60 0.3 404.998 1 5 2.398 0.903 
1088.60 2.0 1088.860 1 20 0.260 
1186.50 2.4 1186.550 1 22 0.050 



Ensemble asteroseismology of ZZ Ceti stars 15 



Table 5 - continued 
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128.10 


3.7 
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1 


3 


0.017 
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1 
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826.40 
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1 


14 
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1 


17 
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1 


19 
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7.55 


288.905 


1 
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0.005 
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206.00 


9.0 
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1 


2 
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0.013 
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2 


7 


1.227 
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268.20 


13.3 


274.612 


1 


3 
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325.50 


13.3 


320.910 


1 


5 


4.590 






353.00 


10.7 


351.912 


2 


11 


1.088 






385.20 


17.0 


382.670 


1 


6 


2.530 






425.80 


13.9 


461.967 


1 


7 


6.167 






533.60 


23.6 


531.504 


1 


9 


2.096 






862.90 


48.9 


862.949 


1 


17 


0.049 




G29-38 


218.70 


1.5 


217.321 


2 


4 


1.379 


2.841 




283.90 


4.8 


282.919 


2 


6 


0.981 






363.50 


4.7 


365.551 
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9 


2.051 
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1150.50 


3.6 


1152.052 


2 


34 


1.552 






1185.60 


3.4 


1185.529 


2 


35 


0.072 






1239.90 


1.9 


1240.220 


2 


37 


0.320 




PG2303+243(®) 


394.4 


7.3 


393.826 


2 


9 


0.574 


0.788 




616.4 


31.4 


616.560 


1 


8 


0.160 






863.8 


7.4 


862.711 


2 


24 


1.089 






965.3 


19.7 


966.590 


1 


15 


1.290 




MCT0145-2211 


462.20 


25 


462.353 


1 


7 


0.153 


1.494 




727.90 


19 


726.912 


1 


13 


0.988 






823.20 


15 


826.663 


1 


15 


3.463 




BPM 30551 


606.80 


11.5 


607.055 


1 


12 


0.255 


0.175 




744.70 


10.5 


744.605 


1 


15 


0.096 




r^T\ no 


1 Qi 1 nn 
ioli.UU 


O.U 


1 Q1 1 nno 


i 


OQ 
ZO 


n nno 
U.UUz 


n nno 
U.UUz 


BPM 24754 


643.70 




643.330 


2 


21 


0.370 


0.938 




1045.10 




1045.204 


1 


20 


0.994 






1234.10 




1234.005 


1 


24 


0.095 






1356.60 




1358.891 


2 


47 


2.291 




KUV 02464+3239 


619.30 


4.0 


618.963 


2 


17 


0.322 


1.640 




777.60 


5.5 


779.541 


1 


12 


1.931 






829.70 


11.6 


829.913 


2 


24 


1.229 






866.20 


9.5 


860.447 


2 


25 


5.704 






993.20 


13.2 


992.707 


1 


16 


0.717 






1250.30 


4.4 


1250.374 


1 


21 


0.121 




PG 1149+058 


1023.50 


10.5 


1023.479 


1 


20 


0.021 


0.021 



BPM 31594 (8) 



401.93 
617.28 



402.453 1 5 0.523 0.321 
617.162 1 10 0.118 
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Table 5 - continued 



Star 


j^obs 


A 


nth 


e 


k 


|A| 








[mma] 








[s] 


[s] 


KUV 11370+4222 


257.20 


5.3 


259.369 


1 


3 


2.169 


0.897 




292.20 


2.5 


291.687 


1 


4 


0.513 






462.90 


3.2 


462.919 


2 


15 


0.019 




HS 1824-6000 


294.30 


8.84 


289.395 


1 


3 


5.005 


2.085 




304.40 


7.66 


301.198 


2 


8 


3.202 






329.60 


13.56 


329.587 


1 


4 


0.013 






384.40 


3.30 


384.520 


2 


11 


0.120 




KUV 08368+4023 


618.00 


16.0 


618.823 


1 


11 


0.823 


0.429 




494.50 


5.5 


494.464 


2 


16 


0.036 




R808 


404.46 


1.99 


400.923 


2 


14 


3.534 


3.499 




511.27 


4.49 


514.497 


1 


10 


3.231 






632.18 


3.41 


629.270 


2 


24 


2.909 






745.12 


3.97 


747.750 


1 


16 


2.630 






796.25 


3.97 


799.402 


2 


31 


3.149 






842.71 


2.81 


844.484 


2 


33 


1.777 






860.23 


3.48 


865.257 


2 


34 


5.030 






875.15 


3.73 


870.376 


1 


19 


4.770 






911.53 


3.19 


913.952 


1 


20 


2.418 






915.80 


5.54 


615.230 


2 


36 


0.573 






952.39 


3.36 


945.909 


1 


21 


6.483 






960.53 


3.68 


967.199 


2 


38 


6.672 






1011.39 


2.54 


1013.941 


2 


40 


2.551 






1040.07 


3.34 


1038.204 


2 


41 


1.866 






1066.73 


2.21 


1066.513 


1 


24 


0.217 






1091.09 


2.36 


1084.277 


2 


43 


6.813 






1143.96 


2.50 


1148.820 


1 


26 


4.860 




G255-2 


685.00 


44 


685.022 


1 


13 


0.225 


0.120 




830.00 


38 


830.218 


1 


16 


0.218 




HLTau-76 


382.47 


16.47 


386.470 


1 


6 


4.001 


2.189 




449.12 


6.7 


447.284 


2 


14 


1.836 






492.12 


7.12 


494.302 


1 


8 


2.182 






540.95 


28.45 


540.790 


2 


18 


0.160 






596.79 


14.40 






10 


1.167 






664.21 


14.94 


663.649 


1 


12 


0.561 






781.00 


9.1 


789.047 


2 


27 


8.047 






799.10 


5.91 


799.328 


1 


15 


0.228 






933.64 


2.40 


933.879 


1 


18 


0.239 






976.64 


6.46 


977.320 


2 


34 


0.680 






1064.91 


11.30 


1064.845 


1 


21 


0.065 






1390.84 


3.92 


1389.281 


1 


28 


1.559 




G232-38 


741.60 


1.9 


741.121 


1 


14 


0.479 


2.155 




984.00 


2.2 


983.679 


1 


19 


0.321 






1147.50 


1.9 


1153.164 


1 


23 


5.664 




G30-20 


1068.00 


13.8 


1068.028 


1 


20 


0.028 


0.028 



(1) Bognar & Paparo (2010), (2) Yeates et al. (2005), (3) Bischoflf-Kim (2009), Vauclair et al. (2000), Voss et al. (2006), 

Pakstiene et al. (2011), Bognar et al (2009), Bradley (1995) 



ment with the spectroscopic value of Tes- To this end, we 
selected the sequence with M* — 0.77Mq and computed an 
additional sequence with Mh = 1.55 x 10~*M«. In this way, 
we obtained a best-fit model with Tes ~ 12 350 K. 

- G185— 32. The pulsation spectrum of this DAV in- 
cludes a period at 215.74 s, quite similar to the dominant 
mode in G117— B15A, but at variance with this star, the dif- 
ference of amplitude between this mode and the remaining 



ones is not so strong in the case of G185— 32. The identifica- 
tion of the £ degree for the periodicities observed in G185— 32 
is not well determined. In particular, the period at ~ 215 s 
is associated with a£=lor£ = 2 mode (Castanheira et 
al. 2004; Yeates et al. 2005). Similarly to G117-B15A, for 
this star the stellar models fit the period at 215.74 s with 
a mode characterized by £ = 1 and k — 2. However, the 
seismological model for this star is more massive than in the 
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Table 6. Structural parameters of the asteroseismological models for the sample of ZZ Ceti stars analyzed in this paper. The quoted 
uncertainties arc the internal errors of our asteroseismic procedure. 



Star logg T.ff M,/Mq Mh/M. Mue/M. Iog(i/L0) log(fl/i?.0) Xc Xq 

[K] 



HS 1531+7436 


8 


28 


± 





06 


12 496 ± 210 





770 


+ 





034 


(1.55 + 


0.23) 


X 


10" 




5 


96 


X 


10" 




-2 


616 


+ 





oil 


-1 


977 


+ 





oil 





332 





655 


GD 244 


7 


97 


± 





04 


12 422 ± 105 





593 


+ 





012 


(1.17 + 


0.36) 


X 


10" 




2 


38 


X 


10" 




-2 


433 


+ 





oil 


-1 


881 


+ 





oil 





283 





704 


G226-29 


8 


28 


± 





06 


12 270 ± 290 





770 


+ 





034 


(2.02 + 


0.31) 


X 


10" 




5 


95 


X 


10" 




-2 


647 


+ 





oil 


-1 


977 


+ 





oil 





332 





655 


HS 0507+0434B 


8 


10 


± 





06 


12 257 ±135 





660 


+ 





023 


(5.68 + 


1.94) 


X 


10" 




1 


21 


X 


10" 




-2 


532 


+ 





021 


-1 


918 


+ 





016 





258 





729 


LP 133-144 


8 


03 


± 





04 


12 210 ± 180 





609 


+ 





012 


(1.10 + 


0.79) 


X 


10" 




2 


45 


X 


10" 




-2 


507 


+ 





010 


-1 


903 


+ 





oil 





264 





723 


EC 11507-1519 


8 


17 


± 





07 


12 178 ± 230 





705 


+ 





033 


(3.59 + 


1.09) 


X 


10" 




7 


63 


X 


10" 




-2 


592 


+ 





021 


-1 


943 


+ 





016 





326 





661 


L19-2 


8 


17 


± 





07 


12 105 ± 360 





705 


+ 





033 


(3.59 + 


1.66) 


X 


10" 




7 


63 


X 


10" 




-2 


602 


+ 





021 


-1 


943 


+ 





016 





326 





661 


GD 66 


8 


01 


± 





04 


12 068 ± 125 





593 


+ 





012 


(4.65 + 


4.37) 


X 


10" 




2 


39 


X 


10" 




-2 


514 


+ 





010 


-1 


896 


+ 





oil 





213 





704 


G132-12 


7 


96 


± 





05 


12 067 ± 180 





570 


+ 





012 


(1.97 + 


0.46) 


X 


10" 




3 


49 


X 


10" 




-2 


486 


+ 





017 


-1 


882 


+ 





014 





301 





606 


G207-9 


8 


40 


± 





07 


12 029 ± 130 





837 


+ 





034 


(4.32 + 


3.50) 


X 


10" 




3 


19 


X 


10" 




-2 


761 


+ 





020 


-2 


017 


+ 





016 





346 





641 


G117-B15A 


8 


00 


± 





09 


11 985 ± 200 





593 


+ 





007 


(1.25 + 


0.70) 


X 


10" 




2 


39 


X 


10" 




-2 


497 


+ 





030 


-1 


882 


+ 





029 





283 





704 


MOT 2148-2911 


8 


05 


± 





04 


11 851 ± 150 





632 


+ 





014 


(7.58 + 


1.79) 


X 


10" 




1 


75 


X 


10" 




-2 


561 


+ 





Oil 


-1 


904 


+ 





Oil 





232 





755 


G38-29 


8 


28 


± 





06 


11 818 ± 50 





770 


+ 





034 


(1.23 + 


0.76) 


X 


10" 




5 


96 


X 


10" 




-2 


716 


+ 





oil 


-1 


979 


+ 





010 





333 





655 


PG 1541+650 


8 


04 


± 





04 


11 761 ± 60 





609 


+ 





012 


(1.56 + 


1.42) 


X 


10" 




2 


46 


X 


10" 




-2 


583 


+ 





010 


-1 


908 


+ 





Oil 





264 





723 


G191-16 


8 


06 


± 





04 


11741 + 90 





632 


+ 





014 


(1.39 + 


0.32) 


X 


10" 




1 


76 


X 


10" 




-2 


590 


+ 





010 


-1 


910 


+ 





oil 





232 





755 


G185-32 


8 


12 


± 





10 


11 721 ± 370 





660 


+ 





023 


(4.46 + 


3.20) 


X 


10" 




1 


22 


X 


10" 




-2 


632 


+ 





051 


-1 


930 


+ 





034 





258 





729 


EC14012-1446 


8 


05 


± 





04 


11 709 ± 95 





632 


+ 





014 


(7.58 + 


2.40) 


X 


10" 




1 


75 


X 


10" 




-2 


583 


+ 





Oil 


-1 


904 


+ 





Oil 





232 





755 


EC23487-2424 


8 


28 


± 





06 


11 700 ± 75 





770 


+ 





034 


(2.02 + 


0.32) 


X 


10" 




5 


95 


X 


10" 




-2 


731 


+ 





010 


-1 


978 


+ 





010 





332 





655 


GD 165 


8 


05 


± 





07 


11635 + 330 





632 


+ 





014 


(7.58 + 


3.28) 


X 


10" 




1 


75 


X 


10" 




-2 


594 


+ 





043 


-1 


904 


+ 





029 





232 





755 


R548 


8 


03 


± 





05 


11 627 ± 390 





609 


+ 





012 


(1.10 + 


0.38) 


X 


10" 




2 


45 


X 


10" 




-2 


594 


+ 





025 


-1 


904 


+ 





015 





264 





723 


HE 1258+0123 


8 


07 


± 





03 


11 582 ± 100 





632 


+ 





014 


(4.46 + 


3.07) 


X 


10" 




1 


76 


X 


10" 




-2 


620 


+ 





014 


-1 


913 


+ 





007 





232 





755 


GD 154 


8 


20 


± 





04 


11 574 ± 30 





705 


+ 





033 


(4.58 + 


1.80) 


X 


10" 




7 


66 


X 


10" 




-2 


705 


+ 





003 


-1 


955 


+ 





003 





326 





661 


GD 385 


8 


07 


± 





03 


11 570 ± 90 





632 


+ 





014 


(4.59 + 


2.86) 


X 


10" 




1 


76 


X 


10" 




-2 


628 


+ 





005 


-1 


962 


+ 





005 





232 





755 


HE 1429-037 


8 


13 


± 





05 


11 535 ± 85 





660 


+ 





023 


(4.68 + 


0.86) 


X 


10" 




1 


22 


X 


10" 




-2 


667 


+ 





018 


-1 


934 


+ 





013 





258 





729 


HS 1249+0426 


8 


02 


± 





02 


11 521 ± 35 





609 


+ 





012 


(3.53 + 


1.08) 


X 


10" 




2 


45 


X 


10" 




-2 


595 


+ 





002 


-1 


896 


+ 





002 





264 





723 


G238-53 


8 


03 


± 





02 


11 497 ± 120 





609 


+ 





012 


(1.54 + 


0.28) 


X 


10" 




2 


46 


X 


10" 




-2 


613 


+ 





002 


-1 


904 


+ 





002 





264 





723 


HS 1625+1231 


8 


02 


± 





04 


11485 + 230 





609 


+ 





012 


(3.52 + 


1.67) 


X 


10" 


-5 


2 


45 


X 


10" 


-2 


-2 


600 


+ 





016 


-1 


896 


+ 





012 





264 





723 


G29-38 


8 


01 


± 





03 


11471 + 60 





593 


+ 





012 


(4.67 + 


2.83) 


X 


10- 


10 


2 


39 


X 


10" 


-2 


-2 


612 


+ 





006 


-1 


901 


+ 





006 





283 





704 


PG2303+242 


7 


88 


± 





07 


11210 + 100 


c 


.525 + 


c 


.12 


(4.54 + 


2.95) 


X 


10" 


-8 


4 


94 


X 


10" 


-2 




.579 + 


c 


.03 


-1 


865 


+ 





032 





279 





709 


MOT 0145-2211 


7 


95 


± 





03 


11439 ± 120 





570 


+ 





012 


(1.43 + 


0.38) 


X 


10" 




3 


50 


X 


10" 


-2 


-2 


573 


+ 





014 


-1 


879 


+ 





012 





301 





686 


BPM 30551 


8 


19 


± 





05 


11435 + 40 





705 


+ 





033 


(4.36 + 


0.26) 


X 


10" 


-6 


7 


66 


X 


10" 


-3 


-2 


714 


+ 





006 


-1 


949 


+ 





006 





326 





661 


GD 99 


8 


01 


± 





13 


11 395 + 25 





660 


+ 





023 


(1.36 + 


0.52) 


X 


10" 


-5 


1 


22 


X 


10" 


-2 


-2 


671 


+ 





005 


-1 


950 


+ 





068 





258 





729 


BPM 24754 


8 


03 


± 





03 


11 390 + 50 





609 


+ 





012 


(4.51 + 


2.72) 


X 


10" 


-6 


2 


46 


X 


10" 


-2 


-2 


626 


+ 





Oil 


-1 


902 


+ 





001 





264 





723 


KUV 02464+3239 


7 


93 


± 





03 


11 360 + 40 





548 


+ 





014 


(4.71 + 


2.45) 


X 


10" 


-8 


4 


21 


X 


10" 


-2 


-2 


579 


+ 





006 


-1 


876 


+ 





006 





290 





697 


PG 1149+058 


7 


94 


± 





02 


11 336 + 20 





570 


+ 





012 


(5.29 + 


2.45) 


X 


10" 


-5 


3 


69 


X 


10" 


-2 


-2 


579 


+ 





001 


-1 


875 


+ 





002 





301 





686 


BPM 31594 


7 


86 


± 





03 


11 250 + 70 





525 


+ 





012 


(5.36 + 


1.87) 


X 


10" 


-5 


4 


93 


X 


10" 


-2 


-2 


545 


+ 





009 


-1 


851 


+ 





009 





279 





709 


KUV 11370+4222 


8 


06 


± 





03 


11 237 + 80 





632 


+ 





014 


(1.40 + 


0.64) 


X 


10" 


-5 


1 


76 


X 


10" 


-2 


-2 


668 


+ 





007 


-1 


911 


+ 





008 





232 





755 


HS 1824-6000 


7 


95 


± 





08 


11234 + 400 





570 


+ 





012 


(1.43 + 


0.62) 


X 


10" 


-5 


3 


50 


X 


10" 


-2 


-2 


605 


+ 





050 


-1 


879 


+ 





030 





301 





686 


KUV 08368+4026 


8 


02 


± 





03 


11 230 + 95 





609 


+ 





012 


(1.42 + 


0.52) 


X 


10" 




2 


45 


X 


10" 


-2 


-2 


646 


+ 





010 


-1 


899 


+ 





007 





264 





723 


R808 


8 


18 


± 





05 


11213 + 130 





705 


+ 





033 


(3.59 + 


1.70) 


X 


10" 




7 


63 


X 


10" 


-3 


-2 


738 


+ 





008 


-1 


944 


+ 





008 





326 





661 


G255-2 


8 


11 


± 





04 


11 185 + 30 





660 


+ 





023 


(4.45 + 


2.12) 


X 


10" 


-6 


1 


22 


X 


10" 


-2 


-2 


709 


+ 





002 


-1 


928 


+ 





003 





258 





729 


HLTau-76 


7 


89 


± 





03 


11111 + 50 





548 


+ 





012 


(1.83 + 


1.03) 


X 


10" 


-4 


4 


19 


X 


10" 


-2 


-2 


579 


+ 





005 


-1 


857 


+ 





005 





323 





697 


G232-38 


7 


99 


± 





04 


10952 ± 120 





593 


+ 





012 


(5.19 + 


1.87) 


X 


10" 




2 


38 


X 


10" 


-2 


-2 


666 


+ 





015 


-1 


888 


+ 





010 





283 





704 


G30-20 


7 


91 


± 





02 


10 950 + 15 





548 


+ 





012 


(5.34 + 


2.18) 


X 


10" 


-5 


4 


20 


X 


10" 


-2 


-2 


618 


+ 





002 


-1 


863 


+ 





002 





290 





697 



case of G117— B15A. For G185— 32 we adopted an astero- 
seismological model that closely fit the period at ~ 215 and 
at the same time it matches the set of observed periods with 
mostly £ = 1 modes. 

- GD 154. This star shows three pulsation modes. The 
mode with period at 402.6 s is a unstable and low amplitude 
mode, as compared with the remaining two modes (PfeifTer 
et al. 1996). Since the amplitude of the long period modes 
(1088.6 and 1186.5 s) are very similar, and since the period 
at 1186.5 s is probably a dipole mode (Pfeiffer et al. 1996) we 
favor models that fit these periods with £ = 1 modes. Gen- 
erally, the solutions have a stellar mass between 0.6323 and 
O.7O5M0 with thin H envelopes (Mh ~ 10"* - 10"^° M,). 
Among them, we choose the solution with M* — 0.705Mq 
and Mh ~ 4.58 x 10~^''M« because it has surface parame- 
ters in agreement with spectroscopy. Other similar solutions 



have Toff ~ 11 200 K, but in these cases the period at 402.6 
s is identified with £ — 2. 

- G238-53, G132-12 and LP 133-144. These three 
stars also have a period near 215 s. In all the cases, this 
mode has an identification k = 2 when £ = 1. This £ and 
k identification is an intrinsic property shared by all the 
asteroseismological models of this study. This can be seen 
from Fig. [S] that shows that the periods with £ = 1 and 
k = 1 are always too short to match the ~ 215 s period. 

- R548. This star has a slightly higher effective tem- 
perature and a spectroscopic stellar mass a bit larger than 
G117— B15A. Frequently, both stars are analyzed together 
due to these similarities and several periods in common. In 
the first attempts to fit the periods of R548 we obtained so- 
lutions with high mass, but were discarded because the 212 s 
periods was identified with £ — 2 according to those models. 
Also, intermediate mass solutions were obtained. Generally, 
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the modes with periods at 318.07 s and 333.64 s are the 
most poorly matched by the models, and they are identified 
with £ = 2. In order to found a best-fit model for this star, 
we were forced to employ several restrictions. We assumed 
that the mode with the period at 212.95 s has 1=1 and 
fc = 2, and fixed 1 = 1 also for the mode with the period at 
274.272 s (Yeates et al. 2005). We found an asteroseismolog- 
ical model with M, = O.6O9M0, larger than the stellar mass 
obtained for G117-B15A (M. = Q.593M0) and a T^s lower, 
in contrast to the trend indicated by spectroscopy and by 
the previous studies (Bradley 1998; Castanheira & Kepler 
2009). However, the surface parameters characterizing the 
best fit model are within the uncertainties of spectroscopy. 

- MCT 2148-2911, PG 1541+650, HE 1429-037 
and HS 1824—600. These are low-mass white dwarfs, with 
spectroscopic masses of 0.515, 0.502, 0.514 and 0.427Mq, re- 
spectively. These values are obtained by extrapolation from 
our evolutionary model grid. However, our asteroseismolog- 
ical models for these stars do not have the lowest mass of 
our model grid (Q.525M0), but instead, they result in inter- 
mediate masses: 0.632, 0.609, 0.660 y O.57OM0, respectively. 
The Toff values of these models are in agreement with the 
spectroscopic inferences. 

- GD 244. For this star we have not been able to found 
any plausible seismological model with an effective temper- 
ature close to the spectroscopic value (Toff = 11680 K). In 
order to adopt a seismological model, we considered that the 
large amplitude modes are 1=1, and found an acceptable 
solution with a Tcff = 12 422 K, markedly higher than the 
spectroscopic one. On the other hand, the gravity and stel- 
lar mass of the adopted seismological model are compatible 
(within the uncertainties) with the spectroscopic estimates. 

- G207— 9. For this DAV we obtain a massive seismo- 
logical solution, with Af» = O.837M0. However, a second 
solution, although with a slightly worse period match, ac- 
cording to $ = 1.496 s, is obtained for a lower mass 
(M» = O.6O9Af0), characterized by a thick H envelope 
(Mh = 1.41 X 10~^M,). A degeneracy of solutions for this 
star has been also found by Castanheira & Kepler (2009). 

- G29— 38. This is a rather pathological case. In spite of 
the fact that this star has T^s ~ 11 800 K, it exhibits a rich 
and complex period spectrum (including 14 genuine eigen- 
modes) which is characteristic of cooler DAVs. Thompson et 
al. (2008), by means of VLT spectroscopy, show that most 
of the periodicities exhibited by this star are 1=1 modes, 
but there are also some 1 = 2 modes and possibly one mode 
with ^ = 3 or 4. However, the seismological solutions for 
this star imply that most of the observed modes should be 
£ = 2. In the asteroseismological model adopted, the only 
£ = 1 mode is associated to the mode with a period 614.4 s 
which has the largest amplitude. 

- PG 2303—243. For this star, the observed modes and 
amplitudes were taken from Pakstiene et al. (2011). These 
authors show that this ZZ Ceti has a very rich pulsation 
spectrum with 24 probably independent modes. However, 
most of these modes show very low amplitudes, below ^ 4 
mma. In our analysis, we considered to be real modes only 
those showing amplitudes higher than ~ 4 mma, leaving 
us with just four periodicities. In particular, we fixed the 
harmonic degree to be ^ = 1 for the two main modes, 616.4 
and 965.3 s, while we allow the remainder modes to be ^ = 1 
or 2. 




M./M^ (spectroscopy) 



Figure 9. Comparison between the value of the stellar mass of 
the 44 DAVs stars analyzed in this work, according to our spectro- 
scopic inference (x— axis) and from our asteroseismological anal- 
ysis (y— axis). The red dashed line represents a perfect match 
between both mass estimates. 



4-3.3 Seismic stellar masses 

In this work, the DA white dwarf evolutionary tracks used 
to derive the spectroscopic masses of the DAVs have been 
employed to infer the asteroseismological masses. Thus, a 
comparison between both sets of values is worth doing. We 
compare in Fig. |§] the spectroscopic and asteroseismologi- 
cal meisses. The dotted line is the 1:1 correspondence. The 
plot reveals that the general agreement between both sets 
of estimations is far from being good, the larger discrep- 
ancies reaching differences up to ~ O.2M0. However, the 
bulk of the points in Fig. |9] accumulate around the dotted 
line, demonstrating that no appreciable offset exists between 
the spectroscopic and asteroseismic estimations of the stellar 
mass. 

The distribution of stellar masses according to astero- 
seismology and spectroscopy is depicted in the histograms 
of the upper and lower panel of Fig. 1101 respectively. The 
mean value of the asteroseismological mass is {Af,)scis = 
0.636 ± O.O19M0, slightly larger (~ 0.95%) than the spec- 
troscopic one, {M.)apcc = 0.630 ±O.O28M0Q. Given the very 
different methods employed to infer both values, the excel- 
lent agreement between these average masses is encouraging. 

Castanheira & Kepler (2008, 2009) have performed the 
first asteroseismological study of an ensemble of ZZ Ceti 
stars. They have studied a total of 83 ZZ Ceti stars including 
the bright variables and also a subset of the SDSS variables. 
The average mass of the ZZ Ceti stars as derived by these 
authors is (M,)sois = O.668M0, about 5% higher than our 
value, (M,)scis = O.636M0. We note that Castanheira & 

^ We do not claim the pulsators are more massive, as there are 
strong selection effects in the search for pulsators. 
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Figure 10. Histograms showing the mass distribution for the 
sample of 44 ZZ Ceti stars considered in this work, according to 
our spectroscopic inferences (lower panel) and our seismological 
analysis (upper panel). 



Kepler (2008, 2009) have included several very massive ZZ 
Ceti stars (M« ^ IMq) that have not been considered in 
our study. Given the fact that the numerical tools used in 
modeling the structure, evolution and pulsations of ZZ Ceti 
stars used by the two groups are independent, and given 
that the samples of stars analysed are not the same, we 
consider that the (M*)sois value derived in this work and 
that derived by Castanheira & Kepler (2008, 2009) are in 
very good agreement. 

4.3.4 The thicknesses of the hydrogen envelope 

One of the most important structural parameters we want 
to constrain through asteroseismology of ZZ Ceti stars is the 
thickness of the H envelope in DA white dwarfs. We have 
found a H layer mass of Mh = (1.25 ± 0.7) x 10"'^M. for 
G117— B15A, about two order of magnitude thinner than the 
value predicted by canonical evolutionary computations, of 
Mh ~ 10~''Mt. Here, the analysis of a large number of ZZ 
Ceti stars allows us to explore the distribution of H enve- 
lope thicknesses from their pulsations. In Fig.[Tl]we present 
histograms of the distribution of H envelope thicknesses. In 
the upper panel we show the results for the complete sam- 
ple of 44 stars. Note that there is a pronounced maximum of 
the distribution for log(MH/M») in the range —5 to —4, al- 
though there exists another, much less notorious maximum 
for log(MH/Mt) between —10 and —9. So, it is apparent 
from the figure that there exists a range of thicknesses of 
the H envelope in the studied DAV stars, with a strong 
peak at thick envelopes and another much lower peak at 
very thin envelopes, and an apparent paucity for intermedi- 
ate thicknesses. In the middle panel of Fig. [TT]we show the 
histogram corresponding to the asteroseismological models 
characterized by canonical (thick) H envelope thicknesses, 
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Figure 11. Upper panel: histogram showing the H envelope 
thickness distribution for the sample of 44 ZZ Ceti stars con- 
sidered in this work. Middle panel: histogram for models with 
canonical (thick) H envelope thicknesses, as predicted by canoni- 
cal evolutionary computations according to the value of the stel- 
lar mass. Lower panel: histogram for models with non-canonical 
(thin) envelope thicknesses, as obtained by means of our artificial 
procedure described in Sect. [2731 



that amount to 11 stars. Finally, in the lower panel we dis- 
play the histogram for the non-canonical thicknesses, that 
is, envelopes thinner than those predicted by standard evo- 
lutionary computations depending on the value of the stellar 
mass. As in previous sections, we refer this kind of envelopes 
as "thin" envelopes. We recall that these "thin" envelopes 
have been generated in this work in order to extend the 
exploration of the parameter space of the models for as- 
teroseismology. Note that in most of the analysed stars (34 
stars from a total of 44) our asteroseismological models have 
"thin" H envelopes, as illustrated in Fig. 1121 It is important 
to note, however, that most of our derived envelope masses, 
even being thinner than the canonical values, cluster close 
to the envelope masses predicted by standard evolutionary 
computations, at variance with those of Castanheria & Ke- 
pler (2009), who found a nearly homogeneous distribution 
of envelope masses in their fits (see their Fig. 8). 

The mean value of the H layer mass is (Mh/M») = 
2.71 X 10~^ according to our results. This value is about 
50 times larger than the value obtained by Castanheira & 
Kepler (2009) with different samples, (Mh/M,) = 5.01 x 
10~^. In spite of this difference, both studies concur to the 
conclusion that an important fraction of DA white dwarfs 
might have been formed with a H mass smaller than the 
value predicted by standard evolutionary computations, a 
conclusion we have already suggested at end of Section 14.21 
on the basis of our results on G117— B15A. 
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Figure 12. The values of the H envelope mass versus the stellar 
mass corresponding to the asteroseismological models of the 44 
ZZ Ceti stars analysed in this work. The thick (orange) curve 
depicts the canonical values for the H envelope thickness. 



5 CONCLUSIONS 

In this paper, we have carried out the first asteroseismolog- 
ical application of the evolutionary DA white-dwarf models 
presented in Ahhaus et al. (2010bj!|. Specifically, we per- 
formed a detailed asteroseismological study of 44 ZZ Ceti 
stars extracted from a sample of bright stars for which the 
surface parameters are accurately known. This sample in- 
cludes the archetypal ZZ Ceti star G117— B15A. The aster- 
oseismological analysis of such a large set of stars has the 
potential to characterize the common properties of the class. 
We have employed a large grid of fully evolutionary mod- 
els characterized by consistent chemical profiles from the 
centre to the surface and covering a wide range of stellar 
masses, thicknesses of the H envelope and effective tempera- 
tures. Our asteroseismological approach represents a signif- 
icant improvement over previous calculations that rely on 
the use of DA white dwarf models characterized by simpli- 
fied chemical profiles at the envelope and/or the core. This 
is the first work aimed at an asteroseismological analysis of 
ZZ Ceti stars that employs fully evolutionary white dwarf 
models. 

Our main results for G117-'B15A are: 

• We found an asteroseismological model for G117— B15A 
with Teft = 11985 ± 200 K, log 3 = 8.00 ± 0.09 and 
M* = 0.593 ± O.OO7M0, in excellent agreement with the 
spectroscopic determinations. 

• For the first time, we break the degeneracy of the as- 
teroseismological solutions for this star reported by previous 

^ Detailed tabulations of the chemical profiles for diff'erent stellar 
masses and effective temperatures are available at our web site 
|http : / / www . f caglp . unlp . edu . ar /evolgroup | 



studies regarding the thickness of the H envelope, depend- 
ing on the fc-identification of the three periods exhibited by 
G117— B15A, although it is fair to say that we are matching 
3 periods by varying 3 parameters. We found the identifica- 
tion fc = 2, 3, 4 as the only possible one in the frame of our 
set of pulsation models. 

• Our best-fit model has a H envelope with Mh = (1.25± 
0.7) X 1Q~^M,, about two order of magnitude thinner than 
the value predicted by canonical evolutionary computations, 
of A/h ~ 10~*M» at this stellar mass value. 

• The value of the thickness of the H envelope of our best- 
fit model is in perfect agreement with the predictions of the 
post-LTP scenario proposed by Althaus et al. (2005b) for 
the formation of DA white dwarfs with thin H envelopes. 

• The luminosity of our asteroseismological model allows 
us to infer a seismological parallax of G117— B15A, that is 
substantially larger than its trigonometric parallax. In agree- 
ment with previous works, we argue that the trigonometric 
parallax uncertainty is larger and the seismological deriva- 
tion of the parallax is robust. 

As for the complete sample of 44 ZZ Ceti stars, our 
main results are: 

• We determined the spectroscopic masses of the 44 stars 
analysed using our DA white dwarf evolutionary tracks. 

• The mean value of the asteroseismological mass is 
(M,>seis = 0.636 ± O.O19M0, slightly higher than our mean 
spectroscopic mass, of (M«)spoc = 0.630 ± O.O28A/0. Given 
the completely different approaches employed to derive both 
values, the agreement can be considered as excellent. 

• Our derived value for (M»)scis is in line with the mean 
mass of DA white dwarfs inferred by Tremblay et al. (2011), 
{Mt)oA = Q.612)Mq, and in good agreement with the value 
derived by Falcon et al. (2010), (M*)da = 0.647+0 0" Afg. 

• There exists a range of thicknesses of the H envelope 
in the studied ZZ Ceti stars, in qualitative agreement with 
the results of Castanheira & Kepler (2009). Our distribution 
of H envelope thicknesses is characterised by a strong peak 
at thick envelopes [log(AfH/Af,) ~ —4.5] and another much 
less pronounced peak at very thin envelopes [log(A^H/Af,) ~ 
—9.5], with an evident paucity for intermediate thicknesses. 

• In most of the analysed DAVs (34 stars from a total of 
44), our asteroseismological models have H envelopes thin- 
ner than the values predicted by standard evolutionary com- 
putations for a given stellar mass. However, our envelope 
masses cluster closer to the canonical envelope masses than 
those of Castanheira & Kepler (2009). 

In closing, we note that Tremblay & Bergeron (2008) 
have studied the ratio of He-rich to H-rich white dwarfs in 
terms of Tcff from a model atmosphere analysis of the in- 
frared photometric data from the Two Micron All Sky Sur- 
vey combined with available visual magnitudes. They found 
that this ratio increases gradually from ~ 0.25 for 15 000 > 
Tcff > 10 000 K to about 0.5 for 10 000 > T^a > 8 000 K 
due to convective mixing when the bottom of the H con- 
vection zone reaches the underlying convective He envelope. 
These authors conclude that about 15% of the DA white 
dwarfs should have H envelopes with log(AfH/A/«) between 
— 10 and —8. The asteroseismological results reported in this 
work point to the existence of large fraction of DAV stars 
with H envelopes thinner than canonical values. In particu- 



Ensemble asteroseismology of ZZ Ceti stars 21 



lar, 5 ZZ Ceti stars analyzed have 10"^'' < Mh/M, < 10"**, 
which represents the 11% of the sample of the studied DAV 
stars. This fraction of stars with very thin H envelopes is 
compatible with the results of Tremblay & Bergeron (2008). 

In a detailed asteroseismological analysis of an ensem- 
ble of ZZ Ceti stars, Castanheira & Kepler (2008, 2009) have 
found that the H envelope of these stars could be within the 
range 3 x 10~^° < Mh/M. < 10"", with an average value 
of {Mu/Mt) = 5 X 10"^. In many respects, the results of 
the present study are in excellent agreement with the pre- 
dictions of Castanheira & Kepler (2008, 2009). Our different 
mean value for the H layer mass, (Mh/M,) = 2.71 x 10~^, 
which is about 50 times larger than that found by those au- 
thors, could be due to the fact that our studies are based 
on completely independent sets of DA white dwarf models, 
different pulsational codes, and different samples of stars. 

All these results reinforce the idea that a non-negligible 
fraction of DA white dwarfs with thin H envelopes could 
exist, rendering as a plausible one the scenario proposed 
by Althaus et al. (2005b) for the formation of DA white 
dwarfs with Mh smaller than predicted by the standard the- 
ory. Hopefully, new asteroseismological analysis on a larger 
number of DAV stars, including the ZZ Ceti stars from the 
SDSS, based on fully evolutionary DA white dwarf models 
with realistic chemical profiles like the ones employed in this 
work, will help to place this idea on a firmer basis. 
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